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In  this  report  I attempt  to  fulfill  three  separate  needs  which 
I perceive  to  exist.  First,  and  foremost,  the  field  of  air  cushion 
landing  teclinology  appears  to  me  as  being  one  which  has  not  attracted 
the  attention,  to  any  great  extent,  of  control  system  engineers.  I 
sincerely  hope  that  this  thesis  generates  some  interest  in  the 
control  community  for  the  problems  are  both  challenging  and  stimu- 
lating. Secondly,  this  report  is  written  to  propose  specific  solu- 
tions to  problems  encountered  in  landing  the  Jindivik  aircraft. 

VMle  the  Jindivik  project  was  dropped  during  the  research  phase  of 
this  report,  there  has  been  renewed  interest  in  recent  weeks. 

Finally,  this  thesis  has  been  undertaken  for  my  own  needs:  not  just 
to  meet  degree  requirements,  but  to  finally  use  some  of  the  techniques 
which  have  been  learned  during  my  coursework. 

Many  people  have  contributed  to  this  report  directly  and 
indirectly.  Technical  information  and  computer  resources  have  been 
provided  by  the  Air  Cushion  Group/AFFDL.  Mr.  James  Steiger,  Major 
Jack  Randall  (CAF),  and  Mr.  Dave  Perez,  in  this  organization,  have 
been  most  helpful  and  cooperative.  Several  members  of  the  faculty 
at  AFIT  have  been  extremely  helpful.  Major  Richard  Potter  provided 
excellent  guidance  and  direction  Mi  the  earlier  stages.  Professor 
John  D'Azzo  has  given  much  needed  encouragement  and  many  helpful 
suggestions  in  the  wording  of  the  report.  A very  thorough  reading 
and  criticism  of  the  rough  draft  of  this  report  was  rendered  by 
Dr.  Pete  Maybeck.  Dr.  Maybeck's  comments  were  Invaluable  and  deeply 
appreciated. 
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The  professionalism,  dedication,  and  patience  shown  by  Ms.  Fonda 
Parrott  in  typing  this  thesis  is  sincerely  appreciated  and  will  long 
be  remembered. 

A final  acknowledgement  of  thanks  is  due  my  family  - Sherry  and 
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Abstract 

The  stability  of  an  aircraft  equipped  with  an  air  cushion  recovery 
system  (ACRS)  was  investigated  both  in-flight  and  during  the  landing 
(slideout)  phase  of  flight.  Digital  computer  simulations  were  used 
to  identify  unstable  tendencies  of  a specific  aircraft  (the  Jindivik 
drone).  Control  system  designs  are  proposed  to  eliminate  the  insta- 
bilities . 

A nonlinear,  six  degree  of  freedom,  aerodynamic  model  is  developed 
based  on  available  wind-tunnel  data.  The  model  is  used  in  a computer 
software  package,  EASY,  to  simulate  the  vehicle  in  flight.  Tne 
vehicle  is  shown  to  be  unstable  in  flight  due  to  lateral-directional 
characteristics  created  by  the  installation  of  the  ACRS.  A time 
optimal  control  design  is  developed  using  a coanda-effect  yaw  thruster. 
The  design  is  verified  by  computer  simulations  which  indicate  that 
the  controller  performs  adequately  for  reasonably  small  deviations 
in  sideslip  angle. 

A model  of  the  air  cushion  recovery  system  is  also  developed. 

The  model  is  generated  by  matching  experimental  data  obtained  from 
an  actual  ACRS  with  ll\  spring/damper  units.  The  model  is  verified 
by  computer  simulation,  and  then  slideout  simulations  are  performed. 
These  simulations  reveal  a lateral-directional  instability  in  the 
slideout.  An  optimal  control  design  is  proposed  and  verified,  after 
which  it  is  shown  that  the  stability  of  the  system  is  dramatically 
enhanced  by  the  addition  of  an  optimal  controller. 
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I.  INTRODUCTION 


An  attempt  is  made  here  to  review  the  general  charac- 
teristics of  conventional  landing  gear  and  air  cushion  landing 
devices.  Following  this  review  is  a description  of  some  of 
the  components  that  ruxve  been  used  (or  proposed ) in  the  air 
cushion  technology  field.  A brief  history  of  air  cushion 
research  is  then  presented,  followed  by  a statement  of  the 
problem  areas  approached  in  this  thesis.  The  final  paragraph 
is  an  overview  of  the  remainder  of  the  thesis. 

Some  Prornlnent  Characteristics  of  Conventional  Landing  Gear 

Conventional  landing  systems  have  been  used  for  so  many  years  that, 
perhaps,  some  of  the  undesirable  characteristics  associated  with  these 
systems  are  traits  that  we  have  unconsciously  resigned  ourselves  to 
live  with.  For  instance,  given  a choice  between  constructing  an 
expensive  two  mile  asphalt  runway,  or  a two  mile  grass-covered  landing 
strip,  one  would  likely  pick  the  latter  if  planes  could  safely  land  on 
the  grass.  Thus,  conventional  landing  systems  have  at  least  one  major 
undesirable  characteristic:  they  require  hard  surfaces  for  taxiing, 
taking  off,  and  landing.  Several  other  characteristics  of  conventional 
systems  are  given  below. 

For  takeoffs,  it  is  highly  desirable  that  low  frictional  forces 
exist  in  the  takeoff  direction.  Also,  high  frictional  forces  perpen- 
dicular to  the  takeoff  direction  aid  in  maintaining  a straight  ground 

1 


m 


track.  Low  forward  friction  and  high  lateral  friction  are  obtained 
with  conventional  landing  systems . 

For  landings,  however,  neither  low  rolling  friction  nor  high 
lateral  friction  is  desirable.  Wheel  brakes  are  used  to  increase  con- 
siderably the  forward  friction,  while  skillful  pilot  technique  is  used 
to  circumvent  the  hazard  created  by  high  lateral  friction  in  crosswind 
landings.  Some  aircraft,  notably  the  B-52,  alleviate  the  high  lateral 
friction  problem  by  adjusting  the  angular  orientation  of  the  landing 
gear  to  match  the  angle  between  the  airplane's  heading  on  final  approach 
and  the  runway  centerline  (the  "crab"  angle). 

Conventional  landing  systems  have  the  very  desirable  quality  of 
providing  absorption  of  the  energy  at  inpact.  The  absorption  is 
achieved  primarily  by  the  use  of  dampers  and  the  natural  springiness  of 
the  tires. 

Both  a weight  and  a drag  penalty  are  incurred  by  conventional 
landing  systems.  The  drag  problem  is  frequently  solved  by  making  the 
gear  retractable,  but  this  adds  considerably  to  the  weight  of  the 
system. 

Previously  mentioned  components  include  wheel  brakes,  dampers,  and 
retractable  landing  gear.  Aside  from  adding  weight  to  the  aircraft, 
these  conponents  also  usually  require  hydraulic  actuators.  Hence,  these 
components  require  a larger  capacity  (and  hence,  heavier)  hydraulic 
system,  and  they  also  greatly  increase  the  complexity  of  the  landing 
gear. 

One  characteristic  of  conventional  landing  systems  is  particularly 
desirable:  most  conventional  systems  have  an  adequate  steering  capa- 
bility. Differential  braking,  along  with  nosewheel  steering  on  most 
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aircraft,  allow  precise  control  on  the  ground. 

While  conventional  landing  systems  will  undoubtedly  be  used  for 
many  years  to  come,  the  several  undesirable  traits  listed  above  do 
exist.  A continuing  effort  should  be  made  to  eliminate  these  short- 
comings, and  the  use  of  air  cushion  technology  appears  to  be  the  most 
fruitful  avenue  to  pursue  at  the  present  time  to  achieve  this  objective. 
Below  are  some  of  the  characteristics  of  air  cushion  devices. 


Some  Prominent  Characteristics  of  Air  Cushion  Landing  Systems 

Possibly  the  greatest  potential  advantage  of  air  cushion  landing 
systems  is  the  large  area  over  which  the  weight  is  distributed.  This 
characteristic  allows  landings  on  grass,  sand,  water,  mud,  as  well  as 
asphalt  or  concrete.  Landings  have  already  been  made  on  these  surfaces 
by  the  LA-4  aircraft  (Ref  19).  The  research  program  involving  the  LA-4 
went  further,  however,  and  demonstrated  an  ability  to  taxi  over  various 
surface  irregularities  such  as  ditches,  stumps,  and  plowed  fields. 

For  takeoffs,  it  has  been  noted  that  rolling  friction  is  low  with 
conventional  landing  systems.  Forward  friction  is  also  very  low  with 
air  cushion  systems.  Unlike  conventional  gear,  however,  the  lateral 
friction  is  also  low  and  this  poses  a problem  since  some  steering 
mechanism  must  be  devised  in  order  to  control  the  vehicle  on  the  ground. 
Systems  have  been  designed  to  allow  differential  braking  by  forcing  only 
a portion  of  the  underside  of  the  air  cushion  system  against  the  ground. 
This  method,  known  as  "pillow  braking",  was  used  by  Bell  Aerosystems  on 
the  LA-4  aircraft  (Ref  1^)  with  apparent  success.  Another  method  consi- 
dered for  controlling  the  ground  track  of  air  cushion  vehicles  is 
differential  thrusting  (Ref  12,  Ref  8).  Still  another  device,  one  which 
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is  considered  in  some  detail  in  this  thesis,  is  the  yaw  thruster. 

This  device  holds  considerable  promise,  as  is  later  shorn . 

High  for, vara  friction  is  desired  when  landing  a vehicle,  as  noted 
earlier,  and  effective  methods  of  achieving  this  are  discussed  in  the 
air  cushion  literature.  Pillow  braking  is  one  such  method.  The 
principle  is  the  same  as  discussed  above,  but  areas  on  the  underside  of 
the  air  cushion  system  are  synmetrically  pressed  against  the  ground 
rather  than  differentially.  This  method  was  tested  successfully  on  the 
LA-4  aircraft  referenced  above,  as  well  as  on  the  CC-115  aircraft 
(Ref  8 :l4l).  Another  method  of  recovery  has  been  proposed,  but  has 
not  yet  been  flight  tested.  This  method  uses  a majority  of  the  under- 
side of  the  air  cushion  device  as  a rubber  skid.  A more  detailed 
discussion  of  this  device  is  presented  in  Chapter  4. 

It  was  mentioned  previously  that  dampers  on  conventional  landing 
gear  provide  good  energy  absorption  at  touchdown.  A well  designed  air 
cushion  system  also  provides  good  energy  absorption,  but  heavy  dampers 
are  not  required. 

The  drag  and  weight  of  conventional  landing  systems  have  been 
noted  to  be  undesirable  traits.  Well  designed  air  cushion  systems  offer 
a significant  reduction  in  weight.  The  change  in  drag  as  a result  of  an 
air  cushion  system  varies  considerably  from  one  design  to  another  in  a 
takeoff  or  landing  configuration,  but  can  be  made  very  small  for  cruise 
configurations . 

Since  an  air  cushion  system  does  not  require  dampers,  gear  retrac- 
tion mechanisms  or  wheel  brakes,  the  complexity  may  be  well  less  than 
that  of  conventional  gear.  It  should  be  noted  that  conventional 

4 ! 

-1 


/ 


retractable  landing  systems  have  been  studied  and  improved  upon  for 
many  years.  Yet,  they  are  very  complex  and  require  a multitude  of 


moving  parts.  To  demonstrate  the  simplicity  of  air  cushion  systems,  as 
well  as  to  present  the  common  components  which  are  referred  to  through- 
out the  remainder  of  this  report,  some  definitions  and  terminology  are 
now  given. 

Air  Cushion  Components 

Air  cushion  systems  are  frequently  categorized  by  function.  The 
name  "Air  Cushion  Landing  System  (ACLS)"  is  used  to  denote  an  air 
cushion  system  which  is  designed  for  takeoffs  as  well  as  landings.  The 
name  is  somewhat  misleading,  since  no  mention  of  takeoff  capability  is 
included,  but  nevertheless  this  is  the  nomenclature  which  is  widely 
used  (Ref  17:6).  When  a system  has  been  designed  for  takeoffs  only, 
the  phrase  "Air  Cushion  Takeoff  System  (ACTS)"  is  used.  Similarly, 
when  a system  has  been  designed  to  be  used  only  during  recovery,  the 
phrase  "Air  Cushion  Recovery  System  (ACRS)"  is  used.  While  the  major 
portion  of  this  thesis  is  concerned  only  with  Air  Cushion  Recovery 
Systems,  the  following  few  paragraphs  discuss  all  three  systems  in 
order  to  place  the  ACRS  in  the  proper  perspective. 

A major  component  of  any  air  cushion  system  is  the  inflatable  bag 
called  the  "trunk".  The  trunk  is  composed  of  flexible  material  which 
is  attached  to  the  underside  of  an  aircraft  fuselage  for  the  purpose  of 
containing  a volume  of  air.  The  overall  shape  of  the  trunk  (when 
Inflated)  is  typically  similar  to  an  elongated  torus,  but  the  shape  is 
sometimes  altered  to  achieve  certain  desired  effects.  Figure  1 depicts 
a Jindivlk  drone  aircraft  resting  on  an  ACLS  trunk  (Ref  15).  A cross 
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sectional  view  of  the  trunk  is  depicted  in  Figure  2 along  with  other 
typical  air  cushion  components.  From  this  figure  it  is  seen  that  the 
trunk  serves  as  a reservoir:  it  receives  a supply  of  air  and  simul- 
taneously exhausts  air  through  nozzles  on  its  underside  as  well  as 
sometimes  through  a pressure  relief  valve.  Before  these  three  com- 
ponents are  discussed,  a description  is  first  given  of  what  is  meant  by 
the  term  "air  cushion". 

The  volume  of  air  enclosed  by  the  underside  of  the  fuselage,  the 
inside  surface  of  the  trunk,  and  the  ground  is  generally  referred  to  as 
the  air  cushion.  Air  may  be  supplied  directly  into  this  volume  as 
shown  in  Figure  2,  or,  the  only  flow  into  the  air  cushion  might  come 
from  the  jets  around  the  underside  of  the  trunk.  A given  design  of  an 
air  cushion  system  may  use  either  or  both  methods  of  supplying  the  air, 
but  the  effect  which  is  achieved  is  essentially  the  same  from  either 
method.  The  air  cushion  serves  to  lift  the  vehicle  and  trunk  by  a 
fraction  of  an  inch  above  the  ground,  which  allows  the  vehicle  to  "ride" 
on  air.  The  friction  coefficient  is  reduced  to  a very  small  value. 

The  air  supply  shorn  in  Figure  2 might  be  provided  by  a separate 
engine/compressor  combination  (as  it  was  for  the  LA-4  ACLS).  It  might  also 
come  from  the  propulsion  engine  (as  has  been  designed  for  the  Jindivik 
ACTS). 

The  nozzles  on  the  underside  of  the  trunk  are  usually  placed  all 
around  the  trunk  for  ACLS  and  ACTS  designs.  For  ACRS  desigps,  however, 

4 there  appears  to  be  a good  argument  for  "lubrication" — placing  nozzles 


in  only  the  forward  portion  of  the  trunk.  Such  a placement  allows  the 
forward  area  to  maintain  a low  friction  coefficient.  Two  specific 
advantages  are  derived  from  a reduced  forward-area  friction.  First, 


the  directional  stability  is  enhanced  since  the  resultant  frictional 
force  acts  more  aft  of  the  center  of  gravity.  Second,  the  forward 
section  of  the  trunk  is  acting  to  prevent  the  large  nose-down  pitching 
moments  that  are  present  during  the  landing. 

A pressure  relief  valve  is  frequently  included  in  the  design  of  an 
air  cushion  system.  Its  purpose  is  to  prevent  the  pressure  in  the  trunk 
from  exceeding  some  predetermined  value.  There  are  several  possible 
reasons  for  including  this  component.  Obviously  there  may  be  structural 
limitations  which  might  be  exceeded  with  too  great  a trunk  pressure. 

This  limitation  is  shared  by  the  ACLS,  ACPS  and  ACTS.  The  relief  valve 
might  also  provide  an  improvement  in  stability  by  reducing  fan  stall 
(Ref  5 :15).  Another  reason  for  including  a pressure  relief  valve  in 
the  design  is  to  provide  additional  damping  during  the  landing  inpact 
(Ref  21:7^).  Finally,  the  use  of  a pressure  relief  valve  might  be  the 
primary  means  of  regulating  the  pressure,  particularly  for  a direct- 
flow  type  of  air  supply  (Ref  17:16). 

One  final  component  of  air  cushion  systems  is  the  brake  tread. 
Generally  the  brake  tread  is  made  of  rubber  and  is  simply  used  to  pro- 
vide high  frictional  forces  during  landing.  The  pillow  braking  system 
discussed  previously  is  illustrated  in  Figure  3.  Here  the  brake  tread 
is  placed  only  on  the  underside  of  the  pillows.  This  design  is  very 
attractive  for  an  ACLS,  while  a simpler  approach  is  often  taken  with 
air  cushion  recovery  systems  since  the  design  objective  for  an  ACRS  is 
to  achieve  a high  frictional  force  which,  in  turn,  leads  to  short 
landing  distances.  Therefore,  the  entire  underside  of  an  ACRS  might 
v/ell  be  brake  tread.  Nozzles  are  drilled  through  the  brake  tread 
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where  necessary.  With  this  type  of  recovers'  system,  the  braking  begins 
as  soon  as  the  vehicle  touches  down  and  continues  until  the  vehicle 
comes  to  rest.  Since  the  aircraft  is  sliding  on  the  brake  tread  for 
the  entire  landing,  the  landing  is  often  termed  the  "slideout." . 

Recent  History  of  Air  Cushion  Systems 

Three  major  programs  have  been  initiated  in  attempts  to  demonstrate 
the  feasibility  of  using  air  cushion  technology  to  replace  conventional 
landing  gear  in  some  applications.  The  first  program  was  undertaken  by 
Bell  Aerospace,  Division  of  Textron,  in  the  mid  to  late  1960's.  The 
first  air-cushioned  takeoff  and  landing  was  performed  during  this  pro- 
gram on  August  4,  1967,  on  a converted  Lake  LA-4  aircraft.  Reference  19 
presents  highlights  of  a test  program  which  Bell  Aerosystems  accom- 
plished  under  contract  to  the  Air  Force  Flight  Dynamics  Laboratory. 

' I 

The  program  appears  to  have  been  highly  successful,  but  the  LA-4  air- 
craft was  in  a rather  lightweight  class. 

It  was  decided  that  another  aircraft  would  be  fitted  with  an  ACLS 
to  demonstrate  the  feasibility  of  air  cushion  technology  applied  to 
larger  aircraft.  Hence,  in  the  early  1970' s an  advanced  development 
program  was  initiated  utilizing  the  CC-115  aircraft.  The  CC-115  has  a 
gross  weight  of  41,000  pounds.  Reference  8 presents  a detailed  report 
on  the  objectives,  scope,  and  initial  test  results  of  the  CC-115 
program. 

A third  program  has  been  underway  since  the  early  1970's.  This 
program  has  been  directed  towards  demonstrating  the  feasibility  of  using 
air  cushion  technology  for  launching  and  recovering  high  performance  jet 
aircraft.  The  test  bed  aircraft  selected  was  the  Australian-made 
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Jlndlvik  drone.  The  program  received  a setback  in  the  Sumner  of  1977> 
when  the  government  of  Australia  decided  against  furnishing  the 
Jindivik  aircraft  for  a test  bed. 


Analysis  of  the  Problem 

While  much  effort  has  been  expended  on  the  three  programs  mentioned 
above,  the  vast  majority  of  this  effort  has  been  directed  along  the 
lines  of  improving  the  air  cushion  systems.  Certainly  there  was  (and 
is)  much  to  accomplish  in  this  area,  as  air  cushion  landing  system 
technology  is  still  in  its  infancy.  On  several  occasions,  however,  it 
appears  that  a comprehensive  examination  of  the  stability  and  control 
strategy  was  sorely  needed.  For  instance,  films  of  taxi  tests  with  the 
Jindivik  drone  aircraft  show  poor  directional  stability,  severe  nose 
down  pitching  tendencies,  and  rather  violent  roll  rates.  A taxi 
accident  occurred  during  the  testing  when  lateral  control  was  insuffi- 
cient during  a braking  run  (Ref  15,  Appendix  I). 

At  the  request  of  Flight  Dynamics  Laboratory  personnel,  tins  thesis 
has  been  directed  at  analyzing  the  stability  and  control  of  the  Jindivik 
aircraft  during  landing.  The  research  documented  by  this  report  has 
come  from  a critical  study  of  the  existing  characteristics  of  the 
Jindivik  aircraft  and,  where  needed,  control  strategies  have  been 
proposed  and  analytically  tested.  Three  specific  problems  are  addressed 
in  the  following  four  chapters:  1)  in-flight  lateral-directional 
stability  with  the  trunk  inflated,  2)  lateral-directional  stability 
during  the  slideout,  and  3)  pitch  stabilization  during  the  slideout. 
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Overview 


Modeling  of  the  physical  situation  is  broken  into  two  segments — 
the  aerodynamic  model  and  the  trunk  model.  Chapter  2 is  devoted 
entirely  to  modeling  the  Jindivik  aircraft  in  flight  with  the  trunk 
inflated.  The  third  chapter  makes  use  of  the  aerodynamic  model  to 
demonstrate  an  aerodynamic  lateral-directional  stability  problem 
created  by  the  ACRS,  and  then  gives  a solution  to  this  problem. 

Chapter  4 derives  an  ACRS  model  which  is  then  incorporated  into  the 
aerodynamic  model.  The  fifth  chapter  verifies  the  model,  analyzes  the 
yawing  and  rolling  stability  under  various  conditions,  and  finally 
proposes  control  strategies  to  improve  the  system.  Chapter  6 summa- 
rises the  findings  and  states  some  conclusions  based  on  these  findings. 
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II.  THE  AERODYNAMIC  MODEL 


Aerodynamic  forces  greatly  affect  the  motion  of  an 
aircraft  during  the  touchdown , as  well  as  throughout  a con- 
siderable portion  of  the  landing  ground  run.  A mathematical 
description  of  these  forces  is  derived  in  this  chapter  based 
primarily  on  wind  tunnel  data  presented  in  A.ppendix  A. 
Gravitational  forces , thrust  forces  and  moments , and  aero- 
dynamic forces  and  moments  are  combined  with  kinematic 
relationships , as  well  as  with  information  on  the  existing 
autopilot , to  derive  what  is  referred  to  throughout  this 
report  as  the  "aerodynamic  model". 

Introduction 

An  aerodynamic  model  is  necessary  for  two  major  reasons: 

I)  it  allows  an  analysis  of  in-flight  stability  deficiencies  with  an 
Air  Cushion  Recovery  System  (ACRS)  installed,  and  2)  the  aerodynamic 
model  allows  computation  of  necessary  aerodynamic  forces  and  moments 
during  the  slideout.  The  second  reason  is  obvious,  but  some  expla- 
nation of  the  first  reason  is  in  order. 

Every  test  program  involving  Air  Cushion  Recovery  Systems  has, 
thus  far,  been  conducted  on  an  existing  aircraft  to  which  an  ACRS 
has  been  added.  Hence,  no  thought  was  given  during  the  original 
design  of  the  aircraft  to  the  control  or  stability  with  an  ACRS 
installed.  It  is  conceivable  that  an  aircraft  will  some  day  be 
designed  from  the  beginning  to  land  with  an  ACRS;  but,  until  such 


time,  whenever  an  ACRS  is  used  to  replace  conventional  gear  on  an 
existing  vehicle,  the  stability  of  the  system  must  be  considered. 

An  analysis  of  the  in-flight  aerodynamic  characteristics  of  the 
Jindivik  aircraft  is  presented  in  Chapter  III.  This  analysis  not 
only  demonstrates  some  of  the  detrimental  effects  that  an  ACRS  can 
produce  when  installed  on  an  existing  aircraft,  but  further,  it 
serves  to  confirm  the  first  reason  given  above  for  the  necessity  of 
an  aerodynamic  model. 

Four  general  areas  are  covered  in  the  development  of  the  aero- 
dynamic model  of  the  Jindivik  drone.  The  first  of  these  is  a 
general  description  of  the  Jindivik,  since  this  aircraft  has  some 
unique  characteristics.  Secondly,  the  equations  of  motion  are  pre- 
sented in  both  the  inertial  frame  and  the  body  frame.  The  equations 
of  motion  require  a mathematical  description  of  the  forces  which  act 
upon  the  aircraft;  thus,  the  third  area  is  a development  of  these 
forces.  The  final  area  consists  of  a description  of  the  autopilot 
which  is  presently  installed  on  the  Jindivik. 

General  Description  of  the  Jindivik  Drone 

Tne  Government  Aircraft  Factories  (GAP)  of  Australia  produced 
the  Jindivik  drone  aircraft  in  the  1950's.  Several  excellent  reasons 
for  selecting  the  Jindivik  as  a test-bed  aircraft  are  given  in 
Reference  17,  and  among  the  most  significant  of  these  reasons  are, 

1)  the  engine  inlet  on  this  vehicle  is  above  the  fuselage  and  hence 
ingestion  of  debris  should  not  be  a problem,  and  2)  the  aircraft 
normally  lands  on  a metal  skid  and  thus  little  modification  of  the 
landing  procedure  is  anticipated. 
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Several  physical  characteristics  of  the  Jindivik  are  presented 
in  Table  I.  These  characteristics  are  used  throughout  this  chapter 
in  the  derivation  of  the  model.  Two  entries  in  Table  I are  of 
special  note.  First,  the  approach  speed  of  130  knots  is  considerably 
higher  than  the  approach  speed  of  either  the  LA -4  or  the  CC-iiS 
aircraft.  This  is  significant  since  important  information  will  be 
obtained  in  the  Jindivik  program  concerning  the  recovery  of  high 
speed  ACRS-equipped  aircraft.  Secondly,  the  lateral  control  consists 
of  ailerons  only,  i.e.,  there  is  no  movable  rudder.  This  fact  is 
significant  in  any  control  system  design,  since  whatever  control 
strategies  are  developed  may  well  be  useful  only  for  the  Jindivik. 
Thus,  a clear  distinction  is  made  throughout  this  report  between 
solutions  which  are  specific  to  the  Jindivik  and  solutions  which  are 
considered  to  be  appropriate  for  a more  general  class  of  problems. 

Equations  of  Motion 

Any  of  the  standard  textbooks  on  aircraft  stability  and  control 
can  be  consulted  for  a derivation  of  the  equations  of  motion 
(References  14,  16,  or  4,  for  instance).  Hence,  the  equations  of 
motion  are  stated  here  with  only  some  preliminary  remarks  concerning 
reference  systems,  notation,  and  assumptions.  The  assumptions  are 
introduced  as  necessary  to  simplify  the  problem. 

The  first  assumptions  made  here  are  those  commonly  made  in 
stability  and  control  analyses,  and  are  not  very  restrictive: 

1)  Newton's  laws  of  mechanics  are  assumed  to  hold,  2)  an  earth- fixed 
reference  frame  is  assumed  to  be  an  inertial  frame,  3)  the  earth  is 
assumed  to  be  flat  over  trajectories  of  the  vehicle  which  are  of 
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Table  I 


SOME  PERTINENT  CHARACTERISTICS  OF  THE  JINDIVIK  DRONE  AIRCRAFT 


Parameter 

Value 

Source 

Landing  Weight 

.2620  lbs 

Ref  5 

Approach  Speed 

130  knots  

Ref  5 

Surface  Area,  S 

76  ft2  

Ref  6 

Wing  Span,  b 

19  ft 

Ref  6 

Chord,  c 

. 4 ft 

Ref  6 

zt,  - coordinate  of  the 
thrust  vector,  z^ 

-.375  ft 

Ref  1 

Center  of  gravity 
location 

, 20%  of  MAC 

Ref  6 

Tail  Area,  At 

13  ft2 

Ref  18 

Tail  Aspect  Ratio 

2.89 

Ref  18 

Moment  of  Inertia 
about  Xb,  I xx 

1150  slug-ft2 

Ref  5 

Moment  of  Inertia 
about  yb , Iyy 

1800  slug-ft2 

Ref  5 

Moment  of  Inertia 
about  z\j , Izz 

2840  slug-ft2 

Ref  5 

Longitudinal  Control 

Conventional  Elevator 

Ref  5 

Lateral  Control 

Conventional  Ailerons 

Fixed  Vertical  Stabilizer 

Ref  5 

Propulsion 

Single  Jet  Engine 

Ref  5 
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Interest,  and  4)  the  aircraft  mass  is  assumed  to  be  constant.  These 
assumptions  greatly  simplify  the  equations  of  motion.  A description 
of  the  earth  frame  mentioned  in  the  second  and  third  assumptions  is 
now  in  order. 

The  earth  frame  consists  of  a right-hand,  3-dimensional  Cartesian 
coordinate  system  whose  origin  is  any  arbitrary  point  on  the  surface 
of  the  earth.  Hie  third  (or  ze)  axis  is  perpendicular  to  a plane 
tangent  to  the  earth's  surface  and  points  towards  the  center  of  the 
earth.  The  first  (or  Xg)  axis  lies  in  an  arbitrary  direction  in  the 
tangent  plane  but  is  fixed  with  respect  to  the  earth.  The  second 
(or  ye)  axis  completes  the  orthogonal  system.  It  is  seen  to  be  con- 
venient in  later  chapters  to  'fix  the  Xg  axis  such  that  it  is  parallel 
to  the  runway,  but  for  this  chapter  and  the  next  the  direction 
remains  arbitrary. 

■To  simplify  the  problem  still  further,  it  is  assumed  that  the 
vehicle  is  a rigid  body.  This  is  also  a commonly  made  assumption, 
although  there  are  many  instances  where  this  is  totally  unrealistic. 
Since  gross  vehicular  motion  is  of  interest  in  this  report,  as 
opposed  to  aeroelastic  phenomena,  and  since  the  aerodynamic  loads 
are  relatively  light  throughout  the  approach  and  landing  phase,  the 
neglect  of  structural  deformation  is  considered  reasonable. 

The  equations  which  describe  the  motion  of  a rigid  body  in  the 
(assumed)  inertial  earth  frame  are  as  follows: 

F = m A (1) 


M = — 

dt 


(2) 
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where 


F - vector  force  sura  acting  on  the  body 
m - mass  of  the  body 
A - acceleration  of  the  body 
M - vector  moment  sura  acting  on  the  body 
H - angular  momentum  of  the  body 


While  mathematically  correct,  the  above  equations  are  not 
efficient  computationally,  and  therefore,  the  body- frame  axis  system 
is  introduced  in  order  to  remedy  this  situation.  This  axis  system 
is  defined  in  the  following  manner  for  the  purposes  of  this  report: 
the  first  (or  Xfc,)  axis  extends  from  the  origin,  the  center  of  mass 
of  the  body,  through  the  nose  of  the  aircraft,  parallel  to  the 
fuselage  reference  line.  The  second  (or  y^)  axis  is  perpendicular 
to  the  first  axis  and  passes  through  the  right  wing,  while  the 
third  (z^  axis  completes  the  right-hand,  orthogonal  axis  system 
towards  the  underside  of  the  vehicle.  Figure  4 illustrates  this 
axis  system,  along  with  notation  for  velocities,  forces,  and  moments. 
With  this  notation  the  vector  equations  (1)  and  (2)  can  be  expanded 
to  the  following  set  of  six  scalar  equations  (Ref  16): 


m(U  + QW 

- RV)  = 

IX 

(3) 

m(V  + RU 

- PW)  = 

IY 

(4) 

m(W  + PV 

- QU)  = 

IZ 

(5) 

PIX- 

' RIxz  + 

QR(IZ  - Iy)  - 

RQ^xz  = 

IL 

(6) 

Qly  + 

PR(Ix  - 

Iz)  - R2IXZ  + 

P2lxz  = 

IM 

(7) 

RIz  - 

■ R^XZ  + 

PQ(Iy  — IX) 

QRIxz  = 

IN 

(8) 
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where 


Ix  = moment  of  inertia  about  xb 
Iy  = moment  of  inertia  about  yb 
Iz  = moment  of  inertia  about  zb 
*xz  = product  of  inertia 

The  Ixy  ana  Iyz  products  of  inertia  are  usually  negligible  since 
most  aircraft  are  symmetric  about  the  xb  - zb  plane  (Ref  4 :5).  As 
a first  approximation  Jn  thir  report,  Ixz  is  neglected  also  since 
no  value  for  this  product  of  inertia  could  be  found.  A more  refined 
analysis  should  include  this  term,  although  its  magnitude  is  usually 
much  less  than  that  of  the  moments  of  inertia. 

The  equations  of  motion  and  the  reference  frames  have  been 
presented,  but  an  additional  set  of  equations  are  valuable  when  a 
vector  in  one  frame  is  to  be  represented  in  another  frame.  The 
desired  transition  from  one  frame  to  another  is  usually  accomplished 
by  an  Euler  angle  transformation.  The  Euler  angles  are  typically 
denoted  by  f,  0,  and  <f,  which  represent  the  yaw,  pitch,  and  roll 
angles,  respectively.  The  angles  must  often  be  determined  numeri- 
cally since  they  are  related  to  the  inertial  angular  rates  P,  Q,  and 
R by  three  coupled,  nonlinear  differential  equations.  These  equations 
are 

i = P + Q sin  $ tan  0 + R cos  $ tan  0 (9) 

0 = Q cos  4>  - R sin  $ (10) 

V = (Q  sin  4>  + R cos  $)  sec  0 (11) 

If  values  of  these  angles  are  known  at  a particular  time , then  the 
appropriate  transformation  matrix  for  going  from  the  earth  frame  to 
the  body  frame  (or  vice  versa)  can  be  determined.  Let  C,  represent 
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a transformation  matrix  such  that  if  a vector  expressed  in  body-frame 
coordinates  is  multiplied  by  this  matrix,  the  product  is  a represen- 
tation of  the  same  vector  in  earth-frame  coordinates.  In  terms  of 


I 


the  Euler  angles,  this  matrix  is 

cos  ¥ cos  0 cos  ¥ sin  $ sin  0 
- sin  ¥ cos  $ 
sin  ¥ cos  0 cos  ¥ cos  $ 

+ sin  ¥ sin  $ sin  $ 
-sin  0 cos  0 sin  $ 


Ce  = 

Lb 


sin  ¥ sin  $ 

+ cos  ¥ sin  0 cos  4> 
sin  ¥ sin  0 cos  $ 

- cos  ¥ sin  $ 
cos  0 cos  $ 


(12) 


Conversely,  given  a vector  expressed  in  earth-frame  coordinates,  one 

can  find  a transformation  matrix,  C13,  which  can  transform  those 

e 

coordinates  to  yield  a body-frame  representation  of  the  vector.  It 
can  be  shown  that 


= 5T 


(13) 


Thus  far,  the  equations  of  motion  have  been  developed  along 
with  the  mathematical  relationships  necessary  to  transform  vector 
coordinates  expressed  in  one  frame  to  another  frame.  In  the  next 
section  the  forces  and  moments  which  appear  on  the  right-hand  side 
of  equations  3 through  8 are  developed  for  the  in-flight  configura- 
tion. 


Determination,  of  Forces 

All  forces  discussed  in  this  report  arise  from  one  of  five 
different  sources:  thrust,  gravity,  control,  aerodynamic,  or  ground 
contact.  Ihese  sources  are  identified  by  the  subscripts  t,  g,  c,  a, 
and  f,  respectively.  For  example,  suppose  it  is  desired  to  discuss 
the  force  produced  by  a control  input  in  the  direction,  then  it  is 


21 


represented  as  Yc.  Similarly,  a torque  arising  from  an  aerodynamic 
source  about  the  axis  is  denoted  as  Na.  A discussion  of  control 
forces  produced  by  the  existing  autopilot  is  given  in  the  last 
section  of  this  chapter,  while  the  ground  friction  and  reaction 
forces  are  covered  in  the  fourth  chapter.  The  remainder  of  this 
section  is  devoted  to  gravitational,  thrust,  and  aerodynamic  forces 
and  moments. 

Gravity . In  vector  notation  the  gravitational  force  can  be 
written  as 


G = mge3  (14) 

where  g is  the  local  gravitational  acceleration,  and  ea  is  a unit 
vector  in  the  direction  of  the  third  earth  axis,  ze.  Ey  the  trans- 
formation matrix,  C^,  this  force  can  easily  be  represented  in  body- 
frams  coordinates: 


Z 

L J 


II 

fel 

b 

L J 

e _ pe 
” Lb 


or,  by  use  of  the  definition  of  C®  in  equation  12, 


V 

-mgsin  0 

Yg 

mgcos  0 sin  $ 

mgcos  0 cos  $ 

(15) 


(16) 


This  equation  can  be  used  directly  for  computation  of  the  gravita- 
tional forces  acting  on  the  aircraft. 

Thrust.  The  Jindivik  has  a single  jet  engine  which  lies  in 
the  Xj-,  - Zb  plane.  The  thrust  vector  is  assumed  to  be  parallel  to 
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the  fuselage  reference  line  but  displaced  from  the  center  of  gravity 
by  a distance  zp.  Urns,  there  is  only  one  force  component  and  one 
moment  component: 


Xp  = Pp 

(17) 

Mp  = Zp  Fp 

(18) 

where  Fp  represents  the  magnitude  of  the  thrust. 

Aerodynamic  Forces  and  Moments.  In  contrast  with  the  straight- 
forward determination  of  the  gravitational  and  thrust  forces,  the 
aerodynamic  forces  are  difficult  to  obtain  with  great  precision. 

Two  primary  sources  are  relied  upon  for  data  which  is  necessary  for 
the  computation  of  these  forces  (Refs  1 , 15 ) , but  the  data  is  less 
than  exhaustive . In  several  instances  extrapolations  are  necessary, 
but  no  claim  is  made  as  to  the  accuracy  of  such  extrapolations. 

A method  commonly  used  to  obtain  expressions  for  the  aero- 
dynamic forces  and  moments  is  that  of  forming  a first-order  Taylor 
series  expansion  about  some  operating  point  (Refs  1^,  16).  This 
method  is  explained  below  in  general  terms,  after  which  an  example 
is  presented. 

Let  F represent  a function  which  is  analytic  in  some  region  of 

interest,  and  furthermore  let  F be  dependent  upon  a set  of  m 

T 

variables:  (A1}  A2,  ...,  Am)  = A.  Define  A*  to  be  the  values  of 
the  dependent  variables  which  make  up  A at  some  point  within  the 


region  of  interest,  and  call  this  point  the  equilibrium  point.  Then 
the  Taylor  series  expansion  of  F to  first  order,  denoted  by  F'(A), 


is  given  by  the  formula 

F(A)  = F'(A)  = F(A»)  + (Ax  - 


A?)H^)  + (Az  _ A*). 


11  A* 


*N3F(A) 

2)uT~ 


+ ...  + (Am  - Aj-jj) 


*s8F(A) 

m;3Am 
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The  wind  tunnel  data  presented  in  Appendix  A is  used  to  provide  the 
value  of  each  function  at  the  equilibrium  points  (i.e.,  the  F(A*) 
values),  while  a combination  of  given,  derived,  and  wind-tunnel  data 
is  used  in  determining  the  partial  derivatives  indicated  in  equation 
19.  Furthermore , while  many  dependent  variables  affect  the  indi- 
vidual forces  or  moments,  physical  insight  and  experience  have  been 
used  in  the  past  to  limit  the  number  of  these  variables  to  those  few 
which  are  most  significant  (Refs  4,  14,  and  16).  A tabular  listing 
of  the  relative  importance  of  different  dependent  variables  can  be 
seen  indirectly  in  Reference  16,  page  4.114,  where  each  of  24 
"stability  derivatives"  are  given  a score  based  on  their  relative 
importance.  Stability  derivatives  are  directly  related  to  the 
partial  derivatives  indicated  in  equation  19.  In  order  to  illustrate 
the  Taylor  series  expansion  method,  an  example  of  a force  expansion 
is  now  given. 

Consider  the  aerodynamic  force  in  the  z^,  direction,  Za.  The 
dependent  variables  which  have  been  determined  to  be  the  most  signi- 
ficant contributors  to  Za  are,  1)  rate  of  change  of  angle  of  attack, 
a,  2)  pitching  velocity,  q,  and  3)  the  elevator  deflection  angle,  <5g. 
Thus,  in  equation  19  the  variables  represented  by  A are 

A = (a,  q,  <5E)  (20) 

The  reader  may  well  wonder  why  angle  of  attack,  a,  and  forward 
velocity,  U,  have  been  omitted  from  the  listing  of  independent 
variables  since  Za  obviously  changes  when  either  of  these  variables 
are  changed.  The  reason  for  the  omission  is  that  these  variables 
are  included  in  the  computation  of  the  equilibrium  forces  in  each 
iteration  of  a digital  computer  program.  In  effect,  a new 


equilibrium  point  is  established  during  each  iteration  of  the 
program,  with  the  Taylor  series  expansion  being  used  only  to  com- 
pute the  effects  of  the  control  forces  and  the  "dynamic  variables". 
Ihe  dynamic  variables  are  defined  in  this  report  as  variables  which 
are  either  inertial  angular  velocities  (P,  Q,  R),  or  rates  of  change 
of  the  attitude  of  the  aircraft  with  respect  to  the  relative  wind 
(i.e.,  rate  of  change  of  angle  of  attack,  a,  or  rate  of  change  of 
sideslip  angle,  g) . 

The  independent  variables  identified  by  equation  20  are  now 
substituted  for  A in  equation  19,  and  F in  equation  19  is  replaced 
by  Za,  to  yield 

3Z 

Za(&>  q,  «E)  = za<*)  + <&  " 

„ n (2D 

3 Zp  w 3 Zp 

+ (q  - q*)g^—  + («E  " 6I>§fc 

d<i  (#)  d0E  (*) 

In  this  equation,  as  well  as  throughout  this  chapter,  the  symbol 
is  used  to  indicate  that  the  function  is  evaluated  at  the  equili- 
brium point.  Also,  throughout  this  report  the  control  variables 
and  the  dynamic  variables  are  assigned  a value  of  zero  at  every 
equilibrium  point.  Hence,  in  Equation  21,  a*  = q*  = <5^  = 0,  and  the 
equation  can  be  written  in  a more  compact  form 

za  = ZaQ  + Zct&  + Zq1  + Z<$e6E  (22) 

where  the  definitions  of  Z^,  Zj,  Zq,  and  Z$E  are  obvious.  If  the 
coefficients  of  a,  q,  and-  6E  are  divided  by  the  moment  of  inertia 
about  the  Z5  axis,  I2Z,  the  familiar  "dimensional  stability  deriva- 
tives" are  obtained  (Refs  4,  16).  A brief  description  of  how  the 
equilibrium  force  Z^  and  the  values  of  Z* , Zq,  and  ZgE  are  actually 
computed  is  presented  next. 
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In  static  wind  tunnel  testing,  as  was  performed  in  Reference  15, 
the  data  is  obtained  by  subjecting  a model  to  the  flow  in  the  tunnel 
and  measuring  the  resultant  forces  while  the  model  is  placed  system- 
atically in  different  attitudes  and  configurations . Ihe  resulting 
data  is  expressed  in  terms  of  "static"  aerodynamic  coefficients 
which  are  defined  as  follows: 


(23) 

CD  = §3 

(24) 

(25) 

vP 

II 

■&I* 

(26) 

o 

r— • 

H 

Itr1 

(27) 

r - N_ 

Cn  qSb 

(28) 

where, 

pVT2 

q = dynamic  pressure,  — — 

S = planform  surface  area 
c = mean  aerodynamic  chord 
b = wing  span 
L = lift 
D = drag 

M = pitching  moment 
Y = side  force 
= rolling  moment 
N = yawing  moment 

It  has  been  necessary  here  to  change  the  symbol  for  rolling  moment 
given  in  Figure  4 (L)  to  Ll}  so  as  to  distinguish  between  the  usual 
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symbol  for  lift  and  the  symbol  for  roll!:. ; moment.  Both  lift  and 
rolling  moment  are  symbolised  by  1 throughout  the  remainder  of  this 
report  whenever  there  is  no  apt  ' jilt  . D . Lrsl  two  coefficients 


above  are  necessary  to  cc  ::pu*.  . 

The  relationship  betwe-  •: . I-  , t.  : 
below.  Clearly,  can  be  computed  as  follows: 

S' i_  “ ^ COS  3 *■*  — n d 

= -qSCCr,  ccs  a + Cp  sin  i) 


..  ”1  — • r-»o  r* 


(29) 


Thus,  a knowledge  of  Cp,  Cp,  q,  and  a is  all  that  la  ne jessary  to 
compute  2^.  Appendix  A presents  analytic  approximations  of  Cl  and 
Cp  (as  well  as  the  other  static  aeradynaric  coefficients ) as  functions 
of  a,  g,  and  the  flap  setting.  Thus , a computer  program  calculation 
of  (or  any  other  equilibrium  force  or  moment)  during  ary  one 
iteration  requires  only  the  flap  setting,  the  true  velocity,  and  the 


orientation  of  the  aircraft  to  the  flow  (l.e.,  angle  of  attack  and 
sideslip  angle).  The  computer  program  used  for  all  simulations  in 
this  report  performs  calculations  of  equilibrium  forces  and  moments 
in  the  manner  just  described. 

Wind  tunnel  data  for  the  dynamic  variable  coefficients  (Z£  and 
Zq  in  Equation  22)  are  not  available.  Approximate  values  for  these 
two  coefficients  are  derived  in  the  second  section  of  Appendix  A. 

Ihe  elevator  control  coefficient  is  also  needed  in  Equation  22, 
and  there  is  sufficient  wind  tunnel  data  to  establish  this  coeffi- 
cient. Both  the  data  (Tables  XI-A,  XII-A,  and  XEII-A)  and  analytic 
approximations  of  the  data  are  given  in  Appendix  A. 

An  example  of  how  one  aerodynamic  force,  Za,  is  conputed  has  been 
given.  Ihe  other  aerodynamic  forces  and  moments  are  conputed  in  the 
same  manner,  but  some  coefficients  in  the  other  expansions  have  been 
taken  directly  from  Reference  1 . These  coefficients  are  given  in 
Thble  XV-A,  Appendix  A.  The  final  step  in  setting  up  the  mathe- 
matical in-flight  model  is  the  conversion  of  the  algebraic  relation- 
ships into  the  conputer  program  which  has  been  previously  mentioned. 
Much  of  the  tedium  of  writing  such  a program  was  avoided  by  the  use  of 
a software  applications  package  created  by  Boeing  Conputer  Services, 
EASY  (Ref  6 ) . A brief  description  of  EASY  and  portions  of  the 
conputer  program  which  have  been  used  to  simulate  and  analyze  the 
Jindivik  are  given  in  Appendix  B. 

Existing  Autopilot . There  are  two  primary  flight  control  sur- 
faces on  the  Jindivik:  the  elevator  and  the  ailerons . Block 
diagrams  of  the  existing  autopilots  which  drive  each  control  surface 
are  presented  in  Figures  6 and  7 • Consider  the  simpler  of  the  two 
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first  — the  pitch  autopilot.  The  pitch  autopilot,  also  referred  to 
as  the  pitch  control  system,  has  the  combined  purposes  of  allowing  a 
ground  crewman  1)  to  command  a desired  pitch  rate  and  pitch  attitude 
and,  2)  to  achieve  and  hold  that  rate  and/or  attitude  until  receiving 
another  command.  The  use  of  pitch  and  pitch-rate  feedback  is  very 
common,  and  furthermore,  it  is  shown  by  simulations  in  Chapter  III 
that  the  pitch  control  system  is  adequate. 

The  roll  autopilot,  or  roll  control  system,  is  similar  in 
structure  to  the  pitch  autopilot.  Roll  angle  and  roll  angle  rate 
can  be  commanded  and  regulated  by  a ground  crewman.  The  outer  loops 
(connected  with  dashed  lines  in  Figure  7)  were  apparently  designed 
to  stabilize  the  heading  of  the  rudderless  vehicle.  Simulations 
were  run  with  the  outer  loops  engaged,  but  this  configuration  was 
found  to  be  'unstable.  Other  simulations  were  made  with  these  loops 
disengaged  (see  Chapter  III),  but  the  vehicle  was  found  to  be 
laterally  unstable  in  this  configuration  as  well.  Since  an  unaccept- 
able situation  is  achieved  in  either  arrangement,  the  outer  loops 
were  eliminated  in  all  subsequent  simulations.  It  is- demonstrated 
in  Chapter  III  that  the  roll  autopilot  using  only  roll  and  roll-rate 
feedback  is  acceptable  if  yawing  motion  is  controlled. 
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III.  DESIGN  AND  ANALYSIS  OP  A YAW  CONTROL  SYSTEM 

The  development  of  an  aerodynamic  model  of  the 
Jindivik  in  Chapter  II  allows  considerable  information  to 
be  obtained  by  performing  dynamic  simulations.  Several 
simulation  results  are  presented  in  this  chapter  which 
serve  to  demonstrate  the  existence  of  a lateral  stability 
problem.  The  instability  is  shown  to  be  one  which  can 
occur  whenever  an  Air  Cushion  Recovery  System  is  installed 
on  an  existing  aircraft.  Since  the  problem  is  of  a 
general  naturet  an  effective  solution  to  the  problem  is 
highly  desirable.  Therefore,  a considerable  portion  of 
this  chapter  is  devoted  to  developing  a yaw  control  system. 
Although  the  system  proposed  in  this  chapter  is  designed 
specifically  to  the  Jindivik , the  basic  structure  can  be 
used  in  a wide  variety  of  situations. 

Introduction 

In-flight  simulations  of  the  Jindivik  aircraft  equipped  with  an 
inflated  Air  Cushion  Recovery  System  are  examined  early  in  this 
chapter  to  identify  any  instabilities,  or  other  undesirable  dynamic 
characteristics,  which  may  be  caused  by  the  installation  of  an  ACRS. 
Simulation  results  are  also  exploited  toward  the  end  of  the  chapter 
to  validate  a yaw  control  system  design.  Since  simulation  results 
are  a vital  part  of  this  chapter,  two  caveats  are  given  at  the  out- 
set. First,  the  scope  of  the  simulations  is  very  limited,  since  the 
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intent  of  this  report  is  not  to  conduct  a detailed  investigation  of 
the  entire  flight  envelope  of  the  Jindivik.  Secondly,  simulation 
results  are  only  as  accurate  as  the  model  which  produces  those 
results.  The  simulations  in  this  chapter  sometimes  reiy  upon  ex- 
trapolated wind  tunnel  data  as  well  as  somewhat  crude  estimates  of 
some  of  the  stability  derivatives.  In  spite  of  these  shortcomings, 
however,  the  simulations  provide  some  very  useful  information,  as 
will  be  seen  shortly. 

The  remainder  of  this  chapter  is  divided  into  three  main  sec- 
tions. The  first  section  consists  of  results  obtained  from  simu- 
lating both  the  bare  airframe  and  the  airframe  and  autopilot 
combination.  The  second  section  develops  a yaw  control  system  to 
stabilize  the  Jindivik  in  flight.  Verification  of  this  yaw  control 
system  is  achieved  in  the  third  section  by  simulating  the  aerodynamic 
model  with  the  yaw  controller  added. 

Simulations  of  Existing  Aircraft  and  Autopilot 

The  first  simulation  results  to  be  considered  are  those  of 
the  bare  airframe.  These  results  serve  as  a baseline  to  which  the 
dynamics  of  all  other  configurations  can  be  compared.  Another 
reason  for  simulating  the  bare  airframe  is  to  isolate  undesirable 
aerodynamic  characteristics  which  may  be  attenuated,  or  even 
exacerbated,  by  the  existing  control  system. 

The  configuration  and  initial  conditions  for  the  bare  airframe 
simulation  are  given  in  Table  II.  A gross  weight  of  2620  lbs  is 
used  since  this  is  a typical  landing  weight.  The  airspeed  is  chosen 
to  be  300  ft/sec  (180  knots),  since  the  present  landing  speed  is 
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Flaps 

ACRS 

Gross  Weight 
True  Velocity,  Vp 
Altitude 
Thrust 

Equilibrium  Pitch  Angle 


Up 

Deployed 
2620  lbs 
300  ft/sec 
Mean  Sea  Level 
807  lbs 
2.07° 


Equilibrium  Elevator  Deflection  3.26° 


217  ft/sec  (130  knots)  and  the  Jindivik  would  fly  a pattern  airspeed 
somewhat  faster.  Equilibrium  conditions  were  computed  by  hand  for 
angle  of  attack,  thrust  setting,  and  elevator  deflection  angle. 
Aileron  deflection  angle  and  rate  are  set  to  zero  for  the  start  of 
the  simulation.  The  flap  setting  is  chosen  arbitrarily  to  be  0°  and 
this  setting  is  maintained  throughout  all  in-flight  simulations.  By 
choosing  an  initial  pitch  angle  equal  to  the  equilibrium  angle  of 
attack  (2.07°),  a velocity  vector  is  established  which  is  parallel 
to  the  earth's  surface.  While  some  lengthy  simulations  were  run 


33 


— 

Simulation  results  for  the  bare  airframe  are  presented  in 
Figures  8 and  9.  The  variables  displayed  in  Figure  8 are  the  Euler 
angles — roll  ($),  pitch  (0),  and  yaw  (y).  These  angles  represent 
the  orientation  of  the  vehicle  'with  respect  to  the  earth.  Roll  angle 
is  the  angle  between  the  yt>  axis  of  the  aircraft  and  the  surface  of 
the  earth,  with  right  wing  down  being  the  positive  direction.  Pitch 
angle  is  the  angle  between  the  vehicle  axis  and  the  surface  of  the 
earth,  with  a nose-up  direction  being  positive.  The  yaw  angle  is  the 
angle  between  the  xb  axis  of  the  vehicle  and  an  arbitrary  Xg  - ze 
plane  in  the  earth  frame,  where  the  positive  direction  for  yaw  is  to 
the  right.*  The  most  notable  aspect  of  these  simulation  results  is 
the  roll  angle — in  less  than  3 seconds  the  roll  angle  reaches  a 
value  of  -100°.  The  pitch  angle  reaches  -6°  in  3 seconds  and  is 
decreasing.  The  yaw  angle  of  the  aircraft  has  veered  -18°  off  the 
original  heading. 

Figure  9 displays  the  time  history  of  the  sideslip  angle,  3, 
and  the  angular  velocity  about  the  z^  axis,  R.  The  sideslip  angle 
is  the  angle  between  the  velocity  vector  of  the  vehicle  and  the 
longitudinal  axis.  The  reason  for  presenting  the  sideslip  angle  is 
that  the  abrupt  changes  in  orientation  of  the  vehicle  shown  in 
Figure  8 are  intimately  related  to  the  sideslip  angle.  The  fol- 
lowing paragraphs  describe  this  relationship  in  some  detail. 

While  most  aircraft  are  symmetric  about  the  xb  - plane,  the 
wind  tunnel  data  presented  in  Appendix  A show  that  this  is  not  true 

*These  descriptions  of  the  Euler  angles  are  given  only  for  the 
benefit  of  the  reader  who  is  unfamiliar  with  flight  dynamics. 

They  are  not  intended  to  be  rigorous  definitions.  Precise 
definitions  can  be  found  in  References  4,  14,  and  16. 
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of  the  Jindivik.  The  side  force  coefficient  (Fig.  A-4),  roll  moment 
coefficient  (Fig.  A -6),  and  yaw  moment  coefficient  (Fig.  A-10),  all 
have  nonzero  values  when  3 = 0.  The  yawing  moment  is  of  particular 
interest  since  the  slope  of  the  yaw  moment  coefficient,  Cn,  versus 
sideslip  angle  is  negative  at  3 = 0°  for  all  values  of  angle  of 
attack.  Whenever  the  slope  of  the  Cn  vs  3 curve  is  negative,  the 
aircraft  is  statically  unstable  (Ref  4 :113). 

Static  stability  can  be  defined  as  a condition  wherein  the 
aircraft  tends  to  return  to  an  equilibrium  point  when  it  has  been 
displaced  from  that  point.  This  definition  is  not  well  suited  for 
nonlinear  systems  since  the  size  of  the  displacement  can  affect  the 
stability.  Thus,  for  the  purposes  of  this  report,  a restriction  on 
the  size  of  the  displacement  is  required.  For  all  values  of  angle 
of  attack,  the  yaw  moment  coefficient  slope  (with  flaps  up)  is  seen 
to  be  negative  within  the  range  -2°<_s<_+2°  (Figs.  A-9  through 
A-12).  To  see  why  a negative  slope  results  in  static  instability, 
it  is  necessary  only  to  consider  the  consequences  of  either  an 
increase  or  a decrease  in  3.  If  3 were  increased  slightly  from  0°, 
a more  negative  yawing  moment  would  be  produced.  But  a more  negative 
yawing  moment  would  eventually  produce  a more  negative  yaw  angle,  and 
consequently  a further  increase  in  3. 

Returning  to  the  simulation  results  in  Figure  9,  the  effect  of 
lateral  static  instability  can  be  seen  directly.  The  values  of 
at  3 = 0°  is  negative.  This  results  in  a negative  torque  about  the 
zb  axis  which  in  turn  produces  a yaw  acceleration  and  rate.  The 
value  of  the  yaw  rate  after  one  second  is  seen  to  be  about  -5*5°/sec. 
A value  of  2°  has  been  reached  by  3 and  consequently  an  even  higher 
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value  of  Cn  is  produced.  After  the  first  two  seconds  the  results  of 
the  simulation  cannot  be  so  easily  interpreted  since  the  relation- 
ships become  very  nonlinear.  However,  some  general  trends  produced 
by  increases  in  sideslip  angle  can  be  noted.  When  S increases,  the 
rolling  moment,  Cl5  becomes  more  negative  as  can  be  seen  in  Figure 
A-o.  The  rolling  moment  produces  a rolling  acceleration  which  brings 
about  a roll  rate  and  a roll  angle  change.  The  abrupt  change  in  roll 
angle  (from  0°  to  -100°  in  less  than  3 seconds)  has  already  been 
noted  in  Figure  8.  The  roll  angle  variation  can  account  for  changes 
in  pitch,  since  the  pitching  moment  (which  acts  about  the  yb  axis  of 
the  vehicle)  depends  upon  angle  of  attack.  A roll  angle  other  than 
0°  results  in  the  lift  vector  of  the  vehicle  no  longer  being  aligned 
with  the  force  of  gravity.  Hence,  an  unbalanced  vertical  force  is 
produced  which  eventually  results  in  a change  in  angle  of  attack. 

The  change  in  yaw  angle  in  Figure  8 is  the  result  of  both  the  yaw 
rate,  R,  and  the  roll  angle,  $.  It  should  be  noted  here  that  all 
variations  in  the  variables  presented  in  Figures  8 and  9 can  be 
explained  in  terms  of  either  the  lack  of  lateral  symmetry  or,  more 
importantly,  the  lack  of  lateral  static  stability. 

Now  that  the  bare  airframe  has  been  simulated  and  discussed, 
simulations  of  the  airframe  augmented  by  the  autopilot  can  be  des- 
cribed. In  order  to  examine  the  pitch  autopilot  and  roll  autopilot 
separately,  the  equations  of  motion  are  decoupled  into  longitudinal 
equations  (equations  involving  forces  along  the  and  axes,  and 
moments  about  y t, ) and  lateral  equations  (forces  along  yt>  and  moments 
about  the  xb  and  z^  axes ) . 


38 


The  pitch  autopilot,  as  mentioned  in  Chapter  II,  makes  use  of 
pitch  angle  and  pitch  angle  rate  feedback.  An  initial  angle  of 
attack  of  3-07°  (1°  higher  than  the  equilibrium  angle  of  attack)  is 
used  to  check  the  longitudinal  dynamics.  This  angle  of  attack  could 
be  considered  to  have  come  about  from  an  instantaneous  gust  along 
the  axis.  A desired  characteristic  of  the  pitch  autopilot  is 
that  it  be  able  to  return  the  aircraft  to  a desired  pitch  angle 
after  a perturbation.  The  purpose  of  this  simulation  is  to  verify 
this  characteristic. 

Figure  10  portrays  the  results  of  this  simulation.  The  initial 
angle  of  attack,  a,  is  seen  to  be  3-07°,  but  in  about  2.5  seconds  a 
has  returned  to  the  equilibrium  value  of  2.07°.  The  pitch  angle 
drops  sharply  and  then  returns  to  a value  of  about  2.04°.  Although 
the  reference  pitch  attitude  is  2.07°,  the  autopilot  has  not  achieved 
this  value  after  5 seconds.  One  possible  explanation  is  that  there 
is  a slight  drop  in  velocity  during  the  simulation,  and  this  changing 
velocity  might  have  produced  the  slightly  low  value  of  pitch  angle. 
The  difference  is  very  slight,  however,  and  the  pitch  autopilot 
appears  to  be  performing  satisfactorily. 

A simulation  of  the  roll  autopilot  is  also  of  interest  as  roll 
control  of  the  vehicle  is  necessary.  Since  the  vehicle  is  statically 
unstable  about  the  yaw  axis,  yaw  rate  and  angle  are  excluded  from  the 
simulation.  This  simulation  is  intended  only  to  verify  that  the  roll 
autopilot  can  restore  and  maintain  a desired  roll  angle  when  all 
states  other  than  roll  angle,  $,  roll  rate,  P,  and  the  roll  auto- 
pilot internal  states  are  fixed. 
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Pig.  10  Verification  of  the  Adequacy  of  the  Pitch  Control  System 


Figure  11  demonstrates  the  adequacy  of  the  roll  axis  autopilot. 
The  autopilot  brings  the  vehicle  to  a 0°  (commanded)  roll  angle  after 
having  been  displaced  by  an  initial  roll  angle  of  5° . The  roll  rate 
is  returned  to  0°/sec  and  the  aileron  deflection  undergoes  a slightly 
underdamped  oscillation  before  settling  down  to  0° . 

A full  six-degree  of  freedom  simulation  is  now  presented.  The 
only  purpose  of  this  simulation  is  to  demonstrate  the  inadequacy  of 
the  existing  autopilot.  The  same  initial  conditions  are  used  for 
this  simulation  as  are  given  in  Table  II. 

The  simulation  results  of  the  airframe  and  autopilot  combination 
are  presented  in  Figures  12  and  13 . Roll  and  pitch  angles  in  Figure 
12  can  be  compared  with  the  bare  airframe  simulation  results  in 
Figure  8.  A much  improved  situation  is  apparently  obtained  in  these 
two  variables,  but  both  angles  are  still  decreasing  after  3 seconds. 
The  peak  yaw  angle  magnitude  is  seen  to  be  even  greater  than  in  the 
bare  airframe  simulation.  Figure  13  reveals  that  the  lateral  static 
instability  exhibited  by  the  bare  airframe  causes  the  deviations . 

Thus,  neither  the  aileron  controller  nor  the  elevator  controller  are 
effective  means  of  controlling  yaw  position  or  rate.  Hie  lateral 
static  stability  of  the  aircraft  has  been  destroyed  by  the  installa- 
tion of  the  ACRS  and  the  existing  autopilot  is  not  equipped  to  handle 
such  a situation. 

Before  proceeding  to  the  next  section,  a list  of  conclusions  may 
prove  useful: 

1.  The  Jindivik  airframe  equipped  with  a particular  ACRS  is 
statically  unstable. 
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2.  The  existing  pitch  controller  adequately  controls  pitch  rate 
and  pitch  angle. 

3.  The  existing  roll  controller  adequately  controls  roll  rate 
and  roll  angle. 

The  performance  of  the  Jindivik  with  the  existing  autopilot 
system  is  inadequate. 

The  use  of  the  term  "adequate"  in  these  conclusions  is  meant  to  imply 
only  that  the  system  does  not  appear  to  be  divergent  for  the  condi- 
tions considered.  Furthermore,  while  a detailed  aerodynamic  analysis 
of  the  Jindivik  is  not  within  the  scope  of  this  report,  a physical 
reason  for  the  statically  unstable  condition  can  be  given. 

■The  vertical  stabilizer  on  conventional  aircraft  provides  ample 
area  aft  of  the  center  of  gravity  to  create  lateral  static  stability. 
An  ACRS  installed  on  the  Jindivik  results  in  a considerable  addition 
of  surface  area  forward  of  the  center  of  gravity  which  offsets,  to 
some  extent,  stabilising  effects  of  the  tail.  Since  this  particular 
stability  problem  can  occur  whenever  an  ACRS  is  placed  on  an  existing 
airframe,  it  is  reasonable  to  seek  a general  solution  to  this  problem. 
While  the  next  section  is  devoted  to  solving  the  lateral  stability 
problem  of  the  Jindivik  aircraft,  the  same  procedures  can  be  applied 
in  many  other  situations. 

Design  of  a Yaw  Control  System 

A fundamental  question  which  must  be  answered  before  beginning 
any  design  process  is,  "What  is  the  purpose  of  the  design?"  In 
light  of  conclusions  stated  in  the  previous  section,  a statement  of 
the  purpose  of  a yaw  control  system  for  the  Jindivik  can  be  made: 


The  purpose  of  the  yaw  control  system  Is  to  provide  lateral  stabili- 
zation of  the  vehicle  such  that  the  vehicle  is  stable  when  the  ACRS 
is  deployed. 

While  many  different  designs  are  possible,  only  three  are  given 
here  as  possible  alternatives.  While  these  three  approaches  are 
certainly  not  exhaustive,  they  serve  to  illustrate  the  wide  range  of 
possibilities.  The  approaches  considered  are  1)  increasing  the  area 
of  the  vertical  stabilizer,  2)  changing  the  aileron  control  law,  and 
3)  devising  a control  law  for  the  yaw  thrusters. 

The  most  straightforward  solution  to  the  problem  is  to  increase 
the  area  of  the  vertical  stabilizer.  The  main  advantages  of  this 
approach  are  simplicity  and  reliability.  The  disadvantages  are 
numerous:  1)  the  added  tail  area  increases  total  drag  throughout 
the  entire  mission,  even  though  the  additional  area  is  needed  only 
briefly  (while  the  ACRS  is  deployed),  2)  structural  modification  of 
the  tail  may  be  rather  expensive,  3)  the  larger  tail  surface  is 
easier  to  detect  by  both  radar  and  visual  means  (with  obvious  impli- 
cations for  military  use),  and  4)  while  the  increase  in  area  in- 
creases stability,  no  means  of  control  is  added. 

Changing  the  aileron  control  law  is  another  alternative  which 
could  be  used  to  meet  the  stated  purpose.  The  approach  would 
essentially  be  that  of  a single  input  (aileron  deflection)  multiple 
output  (sideslip  angle  and  roll  angle)  control  problem.  Modem 
optimal  control  theory  could  possibly  be  applied  to  determine  a 
feedback  law  which  would  stabilize  the  vehicle  and  allow  some 
control  over  both  sideslip  angle  and  roll  angle.  A big  advantage 
is  inexpensive  implementation,  but  there  is  a disadvantage  also. 
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Although,  requirements  of  controllability  are  met  (i.e.,  an  aileron 


deflection  produces  both  a yaw  and  a roll  moment),  the  yawing 
moment  is  only  weakly  affected  by  an  aileron  deflection.  Additional 
analysis  is  required  to  ascertain  the  adequacy  of  this  design 
approach,  however. 

The  third  approach  mentioned  above  is  that  of  designing  a 
control  law  for  yaw  thruster  devices  which  would  stabilize  the 
vehicle.  The  term  "yaw  thruster  device"  has  been  used  to  denote  a 
broad  range  of  mechanisms  which  are  capable  of  producing  a torque 
about  the  zb  axis.  In  order  to  enumerate  advantages  and  disadvan- 
tages, however,  a more  specific  description  must  be  given.  A 
particular  type  of  yaw  thruster  has  been  designed  for  directional 
control  of  the  Jindivik  which  makes  use  of  the  "Coanda  effect". 

The  Coanda  effect  was  discovered  in  1910  by  Henri  Coanda 
(Ref  7 :1).  While  a full  description  of  the  Coanda  effect  yaw 
thruster  can  be  found  in  Reference  7 5 page  1,  a brief  discussion  of 
its  operating  characteristics  is  sufficient  for  the  purposes  of  this 
report . 

Coanda  effect  yaw  thrusters  create  a torque  about  the  vehicle 
Zfo  axis  by  deflecting  a portion  of  the  exhaust  flow  produced  by  the 
jet  propulsion  engine  (see  Figure  14).  The  deflected  portion-of  the 
flow  becomes  attached  to  a metal  plate  on  either  the  port  or  star- 
board side  of  the  exhaust  nozzle.  The  deflection  takes  place  very 
rapidly  and,  therefore,  it  is  assumed  to  occur  instantaneously  for 
the  purposes  of  this  report.  Control  of  the  device  is  achieved  by 
small  solenoid  operated  valves  which  affect  the  pressure  behind  the 
ports  shown  in  Figure  14.  Figure  15  is  a photograph  of  a yaw 
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thruster  installed  on  the  Jindivik  aircraft.  A mathematical  model 


of  the  yaw  thrust  device  which  has  been  designed  for  the  Jindivik 
is  given  later  in  this  chapter,  but  first  the  advantages  and  dis- 
advantages of  this  device  are  presented. 

The  biggest  disadvantage  of  the  coanda-effect  yaw  thruster  is 
that  the  magnitude  of  the  output  of  the  device  is  dependent  upon  the 
thrust  level  of  the  propulsion  engine.  Hence,  for  takeoff  operations 
the  device  produces  its  maximum  torque,  while  during  the  slideout 
only  a fraction  of  the  maximum  torque  is  produced  since  the  thrust 
is  usually  at  an  idle  setting.  During  cruise  portions  of  a mission 
the  thrust  level  will  vary  somewhere  between  these  extremes.  Some 
preliminary  test  data  (Reference  15 ) indicates  that  the  range  of  yaw 
torque  available  is  roughly  100  ft-lb  (at  idle)  to  450  ft- lb  (at 
maximum  thrust). 

Advantages  of  the  yaw  thruster  for  stabilization  of  the 
Jindivik  in  flight  are  numerous.  Since  the  device  is  to  be  used 
for  control  during  the  takeoff  phase,  no  additional  modification  of 
the  Jindivik  is  necessary,  except  for  the  control  circuitry.  It  is 
claimed,  therefore,  that  low  cost  is  an  advantage  of  using  this 
method  for  stabilizing  the  vehicle.  Also,  the  weight  of  added 
control  circuitry  is  much  less  than  the  additional  weight  required 
by  an  increase  in  tail  area.  A final  advantage  is  that  yaw  thrusters 
could  be  used  for  control  of  the  vehicle  as  well  as  for  stabiliza- 
tion. A brief  description  of  how  the  yaw  controllers  could  be  used 
for  control  in  a turn  is  given  later  in  this  chapter. 

Before  any  details  of  the  design  procedure  are  discussed,  a 
mathematical  model  of  the  coanda-effect  yaw  thruster  is  developed. 
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Figure  16  depicts  the  geometry  of  the  situation.  The  yaw  thruster 
produces  a net  force,  Fy^,  parallel  to  the  y-D  axis  by  deflecting 
part  of  the  exhaust  flow  from  the  jet  engine.  The  direction  of  Fyt 
is  positive  when  the  net  force  acts  in  the  positive  y^  direction. 
Hence,  the  following  force  and  moment  equations  can  be  written: 

Yc(yaw  thruster)  = Fyt  (30) 

Nc (yaw  thruster)  = -iyt  Fyt  (31) 

where  Yc  and  Nc  represent  a control  input  force  along  the  y-D  .axis 
and  a control  input  torque  about  the  Z'D  axis  respectively. 


Fyt  ~ yaw  thruster  force  in  y^ 
direction 

lyt  - distance  to  yaw  thruster 
from  vehicle  eg 


Fig.  16  Geometi rj  of  the  Yaw  Thruster  Installed  on  the  Jindivik 


It  should  be  noted  here  that  the  coanda-effect  yaw  thruster  cannot 
provide  a continuous  range  of  forces.  Instead,  it  can  provide  a 
maximum  net  force  in  either  the  positive  yt>  direction  or  the  nega- 
tive yb  direction,  or  no  net  force  in  the  yv,  direction. 
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Since  the  output  of  this  control  device  is  maximum  effort  in  one 
direction  or  another,  it  appears  that  it  should  be  well  suited  for 
use  in  a time  optimal,  bounded  input  control  design.  This  design 
approach  was  undertaken  and  the  main  features  of  this  design  are  now 
presented. 

Analyses  of  the  stability  of  the  Jindivik  early  in  this  chapter 
imply  that  the  negative  versus  B slopes  of  Figures  A-9  through 
A— 13  in  Appendix  A are  responsible  for  unacceptable  performance  in 
the  simulations.  Yaw  moment  coefficient,  Cn,  is  strongly  related  to 
yaw  angle  acceleration,  R.  Furthermore,  sideslip  angle,  B,  is 
closely  related  to  lateral  velocity,  V(6  2 sin-1  (V/VT)),  where  Vt  is 
the  magnitude  of  the  inertial  velocity  of  the  vehicle.  Therefore, 
it  is  reasonable  to  expect  that  the  side  force  and  yaw  moment  equa- 
tions (equations  4 and  8,  Chapter  2)  are  closely  related.  Thus, 
these  two  equations  are  chosen  in  an  attempt  to  simplify  the  problem. 
The  equations  are  repeated  here  for  convenience: 

m(V  + RU  - FW)  = zY  (4) 

RIZ  - PIXZ  + PQ(Iy  - lx)  + QRIxz  = ZN  ' (8) 

It  is  assumed  that  the  pitch  and  roll  autopilots  are  restricting  Q, 

P,  and  P to  small  values  (i.e.,  these  values  are  kept  small  in 
comparison  with  V,  R,  V,  and  R),  such  that  the  following  terms  can 
be  neglected:  WP,  PIxzj  PQ(Iy  - Ix)j  and  QRIxz-  Equations  (4)  and 
(8)  then  become 


m(V  + RU)  = ZY 

(32) 

RIZ  = £N 

(33) 

lb  simplify  the  problem  further,  linear  approximations  for  zY  and 
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IN  can  be  made.  By  techniques  described  in  Chapter  2,  these  summa- 
tions can  be  written  as 

IY  = Ys  + Yt  + Ya  + Yc  + Yf  (34) 

IN  = Nt  + Ma  + Nc  + Nf  (35) 

where  the  subscripts  g,  t,  a,  c,  and  f are  defined  in  Chapter  II. 

For  wings-level,  in-flight  conditions,  equations  34  and  35  reduce  to 

IY  = Ya  + Yc  (36) 

IN  = Na  + Nc  (37) 

The  aerodynamic  side  force  can  be  written  as 

Ya  = qSCy  = qS[CyQ  + Cy-AB]  = qSCyQ  + qSCygAB  (38) 

or  more  compactly, 

Ya  = + YgAB  (39) 

Definitions  of  the  parameters  used  in  equation  38  are 
q = dynamic  pressure 
S = surface  area 
Cy  = side  force  coefficient 

Cy0  = side  force  coefficient  at  an  equilibrium  point 
3Cy 

Cyg  = evaluated  at  an  equilibrium  point 

A similar  treatment  of  the  aerodynamic  yaw  moment,  Na,  yields 

Na  = qSb[Cn(a, S)1  + NraR  + NpAP  + N$aA6a  (40) 

where 


b = 
Cn  = 

Nr  = 

Mp  = 
N«5a  = 


wing  span 

yaw  moment  coefficient 

3^-  evaluated  at  an  equilibrium  point 
3 No 

evaluated  at  an  equilibrium  point 
evaluated  at  an  equilibrium  point 
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Although  Cn  Is  a function  of  both  a and  3,  the  pitch  autopilot  can  be 
assumed  to  keep  variations  in  a small  such  that  Ct^a,  3)  = 0^(3). 
Further  simplification  results  from  a linearization  of  (^(3): 

Cn(3)  = Cj^  + C^^A3  (41) 

where  is  the  value  of  Cn  at  the  equilibrium  point  and  Cn^  is  the 
partial  derivative  of  (^(3)  with  respect  to  3 evaluated  at  the 
equilibrium  point.  Thus,  equation  (40)  can  be  written  as 

Na  = qSbC^  + qSCn.A3  + NrAR  + NpAP  + N6AA<$a  (42) 

Typically  the  last  two  terns  are  small  when  compared  with  the  other 
terns,  so  these  are  neglected.  Then  the  complete  simplified  expres- 
sion for  Na  is  written  compactly  as 

Na  = %o  + NgA3  + wRaR  (43) 

with  the  definitions  of  the  terns  being  obvious  by  comparison  of 
equation  43  with  42. 

The  control  terms  Yc  and  Nc  have  been  defined  in  equations  30 
and  31,  so  the  complete  expressions  for  IY  and  IN  can  be  written: 

SY  = Yao  + YgA3  + Fyt  (44) 

£N  = Nao  + + ~ iytpyt  ' (45) 

Returning  to  equations  32  and  33,  the  side  force  and  yaw  moment 
equations  are  now  approximated  by 

m(tf  + RU)  = Y^  + YgA3  + Fyt  (46) 

RIz  = %<-,  + N3A3  + NrAR  - lytFyt  (47) 

If  the  equilibrium  point  is  specified  to  some  extent,  the  equations 
can  be  further  simplified.  Assuming  that,  1)  U is  nearly  constant, 

2)  the  equilibrium  point  is  always  taken  with  3=0,  and  3)  the 
equilibrium  point  is  always  taken  with  R = 0,  the  following  state- 
ments can  be  made 
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U = Uc,  a constant  (48) 

AB  = 6 = stn-1V/Vip  i (^)V  (49) 

AR  = R (50) 

Rewriting  equations  46  and  47  using  the  stated  equilibrium  conditions 
and  definitions, 

m(V  + RU0)  = + (Yb^)V  + Pyt  (51) 

+ (Ng^)V  + NrR  - lytFyt  (52) 

These  two  equations  can  be  solved  for  the  first  derivatives  of  V and 
R,  producing 


V = (%-)V  + (-Uq)R  + ^Fvf  + — 
mVi>  u m y ^ m 


R = (=%-)V  + &R  + (=i^)Fvt  + 


N. 


kIzVT 


yt 


So 
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In  state  variable  format,  these  equations  are  expressed  as 


V 

an  a12 

V 

t>l 

Cl" 

• 

S5 

R 

+ 

Fyt  + 

R 

a21  a22 

b2 

c2 

(55) 


where  definitions  of  the  terms  in  equation  55  are 


li  l 


iz  1 


mV^ 
NB 
IZVT 


, _ 1 
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a12 Uo 

NR 

w 

_ Ya o 


a22  - j- 


C,  = 


Ziyi 


c2  - j 


m 

Nao 


The  constant  terns  Cj  and  C2  above  represent  the  contributions 
of  Cy  and  Cn  at  3=0°.  These  terms  exist  due  to  either  a lack  of 
symmetry  in  the  one-quarter  scale  model  used  in  determining  the 
wind  tunnel  data,  or  possibly  they  exist  due  to  misalignment  of  the 
model  in  the  wind  tunnel.  Another  possibility  is  that  the  model 
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represents  the  actual  aircraft  very  well,  the  model  was  not  mis- 
aligned, and  the  actual  aircraft  is  really  not  symmetric  about  the 
xb  - Zb  Plane.  As  mentioned  in  Appendix  A,  a 3-inch  pipe  extends 
several  feet  along  the  upper,  port  side  of  the  vehicle  for  the 
purpose  of  supplying  bleed  air  to  the  ACRS.  It  is  conceivable  that 
this  asymmetric  object  produces  the  values  of  CyQ  and  Cj^,  but  it  is 
beyond  the  scope  of  this  report  to  attempt  to  confirm  or  dispute  this 
possibility.  Another  simplifying  assumption  is  now  made:  the 
constant  terms  are  assumed  to  be  zero,  which  is  equivalent  to  saying 
that  the  aerodynamics  of  the  aircraft  are  adjusted  so  as  to  balance 
the  side  force  and  yaw  moment  when  3 = 0. 

With  the  assumption  that  (C1}  C2)  = (0,  0),  equation  55  is  of 
the  form 

x = Ax  + bu  (56) 

where  the  control,  u,  is  bounded.  If  the  problem  is  now  to  find  an 
extremal  control,  u*,  which  transfers  an  arbitrary  initial  state  at 
t = o to  the  final  state  x(tf)  = 0 in  minimum  time,  then  the  problem 
is  referred  to  as  a stationary,  linear  regulator,  minimum-time 
problem  (Ref  13:249)*  Certain  theorems  which  pertain  to  this 
problem  are  proven  in  the  literature  and  presented  here  for  con- 
venience (Ref  13:249). 

Theorem  1 (Uniqueness) 

If  an  extremal  control  exists,  then  it  is  unique. 

Theorem  2 (Number  of  Switchings) 

If  the  order  of  the  system  described  by  equation  56  is  n, 
if  the  eigenvalues  of  A are  all  real,  and  if  a time- 
optimal  control  exists,  then  each  control  component  can 
switch  at  most  (n-1)  times. 


Furthermore,  for  a much  broader  class  of  minimum-time  problems  with 
bounded  controls,  the  form  of  the  control  is  maximum  effort  (bang- 
bang)  throughout  the  interval  of  operation  (Ref  13:2^7).  To  find 
the  optimal  control  solution  to  the  problem  at  hand,  the  optimal 
control  theory  results  stated  above,  and  a linear  transformation  of 
the  state  variables  in  equation  56  are  all  that  are  necessary. 

Although  the  time  response  of  equation  56  can  easily  be  found, 
a convenient  closed-form  solution  for  X2(or  R)  in  terms  of  Xj (or  V)  is 
not  readily  obtained.  Such  a relationship  is  very  useful  for  finding 
an  extremal  control.  A convenient  method  for  obtaining  a closed 
form  solution  is  to  decouple  the  two  states  in  equation  56  such  that 
the  following  form  is  obtained 


1 

►-* 

O 

1 

Zl 

1 

!5 

+ 

z2 

• — 

0 X2 
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where  Xj  and  X2  are  the  eigenvalues  of  A in  equation  56.  The  above 
form  allows  z2  to  be  written  in  terms  of  z;  by  the  following  mani- 
pulations 


Z2  _ X2Z2  + U 

(XjZj  + u^O) 

(58) 

Zl  XjZi  + u 

or 

dz, 

(X2z2  + u)(^-)  = 

(XlZl  + u)g|* 

(59) 

which  implies 

(X2z2  + u)dzi  = 

(XiZ!  + u)dz2 

(60) 

or,  more  conveniently, 

dzj 

dz2 

(61) 

Xizx  + u 

X^Zi  + u 

* 
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Definite  integration  can  be  performed  in  equation  61  from  an  initial 
state,  (zlQ,  z2o)  = zQ: 


1 fZ‘  - 1 f"'2  >A; 

A,*'  A^i  + u a2j  A2C2  + 

Z2„ 


u 


(62) 


'o  o 

Note  that  u is  treated  as  a constant  in  the  above  development.  This 

is  possible  due  to  the  optimal  control  theory  result  stated  earlier 

which  assures  that  the  form  of  the  solution  is  maximum  effort  (i.e., 

constant  along  a trajectory)  throughout  the  interval  of  operation. 

Performing  the  integration  indicated  in  equation  62, 


A lnCAjCj  + u)  1 ln(x2r,2  + u) 

1 zio  2 


(63) 


Z2r 


which,  after  some  algebraic  manipulation,  yields  z2  in  terras  of  zx 
as  desired: 


^2 


= -a-+(s2  + H_)j_a_L^i2Al 

A2  2o  A2  \zi  + u/Ai# 


(64) 


It  is  often  useful  to  examine  the  functional  relationship 
defined  by  equation  64  in  terms  of  the  "phase  plane".  The  phase 
plane  is  a 2-dimensional  representation  of  values  of  z versus  z . 
The  relationship  between  these  two  variables  appears  as  a set  of 
curves  which  are  frequently  referred  to  as  phase  trajectories. 
Reference  10,  pages  470  to  479,  and  reference  9,  pages  447  to  457, 
can  be  consulted  by  the  interested  reader  for  complete  descriptions 
of  linear,  time- invariant  second  order  systems  in-the  phase  plane. 

In  this  report  only  one  of  several  general  cases  is  considered,  that 
being  the  case  of  real  eigenvalues  (A2  and  A2)  with  one  positive  and 
the  other  negative. 
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Fig.  17  Phase  Plane  Trajectories  with  no  Control  Inputs 


Whenever  the  state  variables  of  a second-order  system  are  de- 
( coupled,  as  in  equation  57,  and  when  the  eigenvalues  are  both  real 

with  one  positive  and  one  negative,  the  phase  plane  trajectories  for 
an  unforced  system  have  the  form  depicted  in  Figure  17 . Ihe  origin 


is  referred  to  as  a saddle  point,  and  the  zx,  7.2  axes  separate  the 
phase  plane  into  four  regions  of  interest.  These  regions  are 
Important  because  the  behavior  of  the  system  can  be  determined  from  a 
knowledge  of  the  particular  region  in  which  the  Initial  state  vector 
lies.  For  exarple,  if  the  initial  state  of  the  system  is  such  that  zx 
and  z2  are  both  positive,  the  system  will  diverge  to  the  right  (i.e., 
to  the  condition  z 2 ■*  0 and  zx  ->■“).  It  is  important  to  note  that 
the  trajectories  are  centered  at  the  origin  for  the  unforced  case, 
i.e.,  for  u = o in  equation  6*1.  This  fact  can  be  seen  by  inspection 
of  equation  6h  directly.  Also  of  interest  is  the  shape  of  the  tra- 
( jectories  when  u is  not  zei’o  but  seme  arbitrary  constant. 
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Rearranging  equation  6^  slightly  is  helpful  in  determining  the 
shapes  of  the  forced  trajectories: 

<**  + (fHd&f1) X2/X'  <®  I 

The  form  of  this  equation  suggests  that  an  axis  translation  could 
be  helpful.  Letting  z'2  = (z2  + W>i)  and  z'\  = (z\  + uAi), 
equation  65  becomes 

z2  = (z2q  + u/X2)  (z^-)A2Ai  (66) 

The  form  of  this  equation  is  the  same  as  that  of  6H  except  that  the 
constant  term  multiplying  (s'i)A2^Alis  different.  Thus,  these 

. 

trajectories  have  the  same  overall  trends  as  those  in  Figure  17, 
but  the  exact  shapes  and  the  location  of  the  origin  of  these  tra- 
jectories will  vary  depending  upon  the  value  of  u. 

Consider  now  the  trajectories  in  Figure  18,  where  families  of 

* 

trajectories  are  sketched  for  three  values  of  u:u  = +8,  u = 0,  and 
u = -ft.  The  term  ft  is  used  here  to  denote  the  maximum  magnitude  of 
the  control  u.  Also  shorn  in  Figure  18  is  a portion  of  two  parti- 
cularly important  trajectories  (A  - 0,  B - 0) — those  which  pass 
through  the  origin.  Since  the  optimal  control  solution  must  take 
the  state  of  the  system  to  the  origin,  these  two  curves  constitute 
the  final  trajectory  of  all  solutions.  Furthermore,  since  the  opti- 
mal solution  can  have  at  most  (n  - 1)  switches  of  the  control 
element,  the  only  task  remaining  is  to  determine  a scheme  which  will 
take  an  arbitrary  initial  state  to  one  of  these  two  trajectories  with 
no  control  switching.  When  such  a trajectory  is  found,  the  solution 
is  complete,  for,  as  stated  previously  in  Theorem  1,  if  an  extremal 
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Fig.  13  Phase  Plane  Trajectories  with  Constant  Control  Inputs 


control  exists  it  is  unique.  From  the  sketch  the  strategy  to  be 
followed  is  obvious — for  an  initial  state  to  the  left  of  the  com- 
bined curves  A - 0 - B,  a positive  value  of  u,  +n,  should  be  appli 
while  for  an  initial  state  to  the  right  of  A - 0 - B,  a negative 
value  of  u,  -n,  should  be  applied.  Examples  of  this  strategy  are 
shorn  in  the  figure  by  the  trajectories  which  began  at  Zq  and  z^. 

A mathematical  statement  of  the  foregoing  is  now  developed. 

Let  the  combined  curves  A - 0 and  0 - B be  denoted  by  S(zi  ). 

o 

In  order  to  determine  S as  a function  of  the  states,  consider  'the 
set  of  initial  conditions  such  that  (z\,  z2)  - (z1q,  z2q)  = 0>  and 
u = ±n.  Then  equation  65  becomes 


which  when  solved  for  z2q  yields 


A more  explicit  determination  of  S ( zi  Q ) can  be  made  only  after 
values  of  A j and  A2  are  known . It  is  a simple  matter  to  define  S, 

given  values  for  ft.  A*,  and  A2>  as  is  shown  later.  With  S(z1q) 

determined,  the  optimal  control  strategy  is 

u = +ft  : for  z0  on  A - 0 portion  of  S 

u = +ft  : for  all  z0  such  that  z2q  > S(zj  ) 

u = -ft  : for  z0  on  B - 0 portion  of  S (69) 

u = -ft  : for  all  zQ  such  that  z2q  < S(z1q) 
u = 0 : if  (zio,  z2o)  = (0,  0) 

The  technique  which  has  just  been  developed  is  now  used  to 
develop  an  optimal  control  lav;  for  the  Jindivik.  The  required  steps 
are:  1)  decoupling  of  the  state  equations  (see  equation  55)  into 
the  form  of  equation  57,  and  2)  derivation  of  the  control  law  by  use 
of  equations  68  and  69. 

VJhen  all  values  for  the  elements  in  the  A and  b matrices  of 
equation  56  are  computed,  the  resulting  equation  is 


V 

-.188  -5.23 

v 

.0123 

= 

+ 

R 

-.735  -.267_ 

R 

-.2018 

The  values  of  A were  derived  in  two  ways.  First,  the  values  were 
computed  directly  from  the  definitions  determined  earlier  in  this 
report  for  the  elements  of  A.  These  values  were  checked  by  the  use 
of  the  "LINEAR  ANALYSIS"  feature  of  EASY  (described  in  Appendix  3), 
and  all  values  were  found  to  agree  within  If.  The  b matrix  is  com- 
puted from  the  definitions  for  the  elements  of  the  matrix  which  are 
given  immediately  following  equation  55. 
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Verification  of  the  simplifying  assumptions  made  in  the  deri- 
vation of  this  simple  state  model  is  achieved  to  some  degree  by 
comparing  the  eigenvalues  of  the  full  16-state  model  with  those  of 
the  A matrix  of  the  2-state  model.  Two  eigenvalues  of  the  full- 
state  model  closely  correspond  to  the  2-state  eigenvalues,  as  shown 
below: 


2-State  Model 

Pull-State  Model 

Percent 

Eigenvalues 

Eigenvalues 

Deviation 

Xi 

1.734 

1.677 

3.4/5 

x2 

-2.188 

-2.302 

5.0 % 

The  switching  curve,  S,  is  now  defined  for  this  specific  problem: 


zio  < 0 


Zj  > 0 


The  optimal  control  strategy  is  now  explicitly  determined: 


z2  = 40  '•  u = +9, 


z2  - z2 


z2  > S 


Zo  < S 


u = -fl 


U = +Q 


u = -ft 


Since  the  system  is  unstable  and  the  control  is  bounded,  there 


are  initial  states  from  which  it  is  impossible  to  reach  the  origin. 


The  location  of  these  states  is  depicted  by  the  shaded  areas  in 


Figure  18.  Mathematically,  a solution  to  the  time  optimal  control 


problem  is  seen  to  be  impossible  whenever  the  initial  condition. 


z1q,  and  the  magnitude  of  the  control,  u,  are  such  that  the  values 


inside  the  parentheses  in  equations  71  and  72  are  negative. 


There  are  many  ways  to  implement  the  control  law  stated  in 


equation  7^,  but  for  the  purposes  of  this  report  the  most  direct 


and  straightforward  method  of  implementation  is  sufficient.  Figure 


19  depicts  a block  diagram  of  the  simulation  model  with  the  yaw 


control  system  added. 


The  same  conditions  as  are  used  for  the  bare  airframe  and  air- 


frame with  autopilot  simulations  (see  Table  II)  are  used  for  the 


ya w control  simulations.  The  purpose  of  these  simulations  is  to 


verify  that  the  design  of  the  yaw  controller  is  sufficient  to 


stabilize  the  Jindivik  in  flight.  Since  the  purpose  here  is  to 


demonstrate  the  validity  of  the  design,  only  a limited  number  of 


simulations  are  presented. 


-J 


AIRCRAFT  DYNAMICS 


i 


The  magnitude  of  the  thrust  available  from  the  yaw  thrusters  is 
assumed  to  be  50  lbs.  It  is  acknowledged  that  this  value  may  be 
somewhat  optimistic  based  upon  the  information  given  in  Reference  15, 
but  the  purpose  here  is  merely  to  verify  the  design  procedure.  Addi- 
tional simulations  can  be  performed  after  more  extensive  information 
is  obtained  on  the  Jindivik  yaw  thrusters. 

The  simulations  presented  here  cover  four  different  situations 


as  shown  below' : 

Simulation 

Number 

Initial  Sideslip 
Angle,  e(°) 

Status  of  Lateral- 
Directional  Coefficients 

1 

+2 

Symmetric 

2 

+2 

Unsymnetrlc 

3 

+7 

Symmetric 

4 

+7 

Unsymmetric 

The  meaning  of  the 

terms  ’’symmetric," 

and  "unsymmetric , " is  that  the 
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curves  representing  the  lateral-directional  coefficients  Cy,  Cj,  and 
Cn  are  all  shifted  so  as  to  pass  through  the  origin  in  the  symmetric 
simulation,  while  these  curves  represent  the  actual  wind  tunnel  data 
in  the  unsymmetric  simulations.  Although  the  controller  is  designed 
for  the  symmetric  case,  the  simulations  demonstrate  that  the  con- 
troller functions  adequately  even  when  the  vehicle  is  in  an  unsym- 
metric configuration,  so  long  as  the  initial  sideslip  angle  is  not 
too  large. 

Initial  conditions  of  2°  and  7°  for  sideslip  angle,  6,  are  used 
to  illustrate,  1)  the  effects  of  the  nonlinear  Cn  vs  3 characteristics 
on  the  performance  of  the  yaw  controller,  and  2)  that  there  is  an 
upper  limit  on  the  value  of  sideslip  angle  which  the  controller  can 
tolerate. 

The  first  simulation  results  are  presented  in  Figures  20  through 
22.  Figure  20  demonstrates  the  immediate  return  of  sideslip  angle, 

3,  to  a stable  value  of  0°  and  the  accompanying  build-up  of  yaw  rate. 
This  figure  illustrates  the  dramatic  improvement  of  the  stability  of 
the  system  accomplished  by  the  yaw  controller.  The  basic  airframe 
simulation,  and  the  airframe/autopilot  combined  simulation,  (Figures 
9 and  13)  indicate  a value  of  3 of  7-5°  and  12.5°  respectively  at 
t = 2 sec,  whereas  the  sideslip  angle  with  yaw  control  is  zero  for 
all  t > 1.  sec. 

Figure  21  is  presented  here  to  demonstrate  that  the  orientation 
of  the  aircraft  with  respect  to  the  earth  is  well  behaved.  Unlike 
the  simulation  results  in  Figures  8 and  12,  the  roll,  pitch,  and  yaw 
angles  are  reasonably  constant  with  the  yaw  control  system  engaged. 
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Histories  of  Sideslip  Anple  and  Yai 


ti* 


Figure  22  depicts  the  phase-plane  trajectory  for  this  simulation. 
The  variables  z2  and  z\  in  this  figure  represent  elements  of  the 
vector  z which  is  defined  earlier  in  the  design  phase  by 

z = T^x 

The  initial  state  of  z,  zQ,  is  below  the  switching  curve,  S,  and 
therefore  a negative  yaw  thruster  force  is  commanded.  The  state 
changes  along  a trajectory  towards  the  switching  curve  until  the 
switching  curve  is  encountered,  at  which  time  the  yaw  thruster  force 
is  switched  to  a positive  value.  While  the  2-static  linear  model  is 
designed  to  yield  no  overshoot  of  the  origin,  the  16-state  nonlinear 
model  does  overshoot  the  origin  slightly. 

The  results  of  the  second  simulation  are  presented  in  Figures 
23  through  25.  The  unsymmetric  model  is  used  in  this  simulation  to 
illustrate  the  effects  of  non-zero  values  of  Cy,  Cj , and  Cn  at  3=0°. 
Although  a stable  situation  is  reached  at  t = 2.6  sec,  the  constant 
side  force,  rolling  moment,  and  yawing  torques  significantly  reduce 
the  effectiveness  of  the  yaw  control  system.  The  orientation  of  the 
aircraft,  as  depicted  in  Figure  24,  appears  to  be  entirely  acceptable 
in  this  simulation.  Figure  25  is  the  phase-plane  portrait  for  this 
simulation. 

For  simulation  number  3 a larger  initial  value  of  3 is  chosen 
such  that  the  Cn  versus  3 curve  is  well  out  of  the  linear  range 
(-2°  <_  6 £ 2° ) . The  nonlinearity  does  not  significantly  degrade  the 
effectiveness  of  the  yaw  control  system,  as  3 steadily  decreases  to 
0°  and  both  3 and  R are  stabilized  with  2.5  sec.  The  aircraft 
orientation  can  be  seen  to  be  well  controlled  in  Figure  27.  The 
phase  plane  portrait  is  shown  in  Figure  28. 
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The  final  simulation  results  are  presented  in  Figures  29  through 
31.  An  initial  value  of  8 = 7°  is  used  for  this  simulation,  as  it  is 
for  simulation  number  3,  but  the  unsymmetric  model  is  simulated  here. 


The  shift  in  the  Cn  versus  3 curves  (Figures  A-9  through  A-12)  are 
such  that  a higher  value  of  Cn  is  prevalent  for  any  given  value  of  3. 

This  increase  is  too  large  for  the  yaw  control  system  to  tolerate, 
even  though  the  correct  control  is  commanded  by  the  control  system 
(u  = -n).  Sideslip  angle  begins  to  decrease  in  Figure  29,  but 
upon  reaching  a value  of  approximately  5.5°  the  value  of  Cn  is  near 
its  minimum  value  (see  Figure  A-10).  At  the  value  of  velocity  chosen 
for  this  simulation,  300  ft /sec,  the  minimum  value  of  Cn  is  approxi- 
mately given  by 

Cnmin  " a = 2.2°  * (107)  (76) (19) (-.004)  = -6l8  ft-lb  (76) 

which  is  118  ft-lb  greater  in  magnitude  than  the  maximum  magnitude 
of  the  yaw  controller  (500  ft-lb).  Figures  30  and  31  further  des- 
cribe the  unfortunate  situation  which  results  from  having  too  little 
control  authority  to  meet  the  destabilizing  yaw  characteristics  of 
the  Jindivik. 

Conclusions 

The  simulation  results  depicted  in  Figures  8,  9,  12,  and  13 
vividly  illustrate  a lack  of  stability  in  both  the  bare  airframe  and 
the  airframe  combined  with  the  existing  autopilot.  There  is  certainly 
no  reason  to  have  expected  the  existing  autopilot  to  cope  with  the 
new  aerodynamic  characteristics  which  have  resulted  from  the  addition 
of  an  ACRS.  The  existing  autopilot  performs  well,  however,  if  the 
yaw  is  controlled,  as  illustrated  by  the  simulation  results  in 
Figures  10  and  11. 
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To  control  the  yaw,  a yaw  control  system  has  been  designed 
which,  when  added  to  the  original  autopilot  design  (original  except 
for  the  R,  /R,  and  //R  feedback  loops  discussed  in  Chapter  II) 
stabilize  the  aircraft  at  one  specific  set  of  conditions.  Limita- 
tions do  exist,  however,  as  depicted  in  Figures  29  through  31,  but 
if  the  magnitude  of  6 is  kept  relatively  small  (-2°  <_  8 <_2°)  the 
yaw  controller  performs  well. 

Although  the  simulation  results  are  encouraging,  it  should  be 
noted  tliat  a full  range  of  operating  conditions  must  be  examined 
before  an  implementation  is  attempted.  Particularly,  the  performance 
of  the  yaw  control  system  in  turns  and  at  different  airspeeds  is 
likely  to  be  somewhat  degraded.  Since  there  is  a specific  angular 
velocity,  R,  associated  with  a given  bank  angle  in  a steady  state 
turn,  a provision  for  achieving  a value  of  R should  be  included  in 
an  implementation  design.  Variation  in  air speed  will  probably  be 
less  of  a problem  for  the  yaw  control  system  since  the  variation 
cannot  be  very  much  (from  a typical  traffic  pattern  speed  of  300 
ft/sec  to  a typical  touch-down  speed  of  200  ft/sec). 

Other  applications  of  the  yaw  control  system  design  discussed 
in  this  chapter  are  possible.  While  the  Jindivik  is  a somewhat 
unusual  aircraft  due  to  its  lack  of  a movable  rudder,  aircraft  which 
have  conventional  rudder  systems  might  also  profitably  use  a yaw 
thruster  control  system: . The  addition  of  an  ACRS  could  produce  un- 
desirable aerodynamic  characteristics  which  are  too  severe  for  the 
rudder  alone  to  cope  with.  The  design  procedure  herein  would  need 
to  be  modified  somewhat  to  account  for  the  contribution  of  the 
rudder,  however.  Another  possible  application  of  the  yaw  control 
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design  discussed  herein  is  in  the  case  of  vertical  takeoff  and 
landing  (VTOL)  aircraft  and  short  takeoff  and  landing  (STOL) 
aircraft.  These  types  of  aircraft  need  yaw  control  at  very  low 
airsoeeds  where  conventional  means  of  yaw  control  (rudders)  a^e 
ineffective. 

The  major  contributions  of  this  chapter  are,  1)  the  identifi- 
cation of  an  aerodynamic  problem  caused  by  the  addition  of  an  ACRS, 
and  2)  a proposed  solution  to  this  problem  for  a specific  aircraft. 
Problem  areas  might  also  exist  in  the  slideout  phase  of  operation; 
thus.  Chapters  IV  and  V are  devoted  to,  1)  building  a landing  model, 


and  2)  investigating  the  Jindivik's  stability  during  the  slideout. 


IV.  MODELING  OP  THE  ACRS  TRUNK 

The  aerodynamic  model  developed  in  Chapter  II  pro- 
vides a means  of  computing  aerodynamic  forces  during 
simulations  of  the  slideout.  These  aerodynamic  forces 
constitute  an  appreciable  portion  of  the  total  force  on 
the  vehicle  at  the  beginning  of  the  slideout.  Ground 
reaction  and  frictional  forces , however , predominate 
toward  the  end  of  the  slideout  where  the  dynamic 
pressure  is  small.  Thus,  this  chapter  is  devoted  to 
developing  a mathematical  model  capable  of  computing 
ground  friction  and  reaction  forces  which  act  through 
the  trunk  of  an  ACRS. 

Introduction 

Several  different  approaches  to  developing  a trunk  model  have 
been  used  in  the  past.  Ihe  simplest  models  have  assumed  the  trunk 
to  act  as  a set  of  two  linear  torsional  springs  (one  for  roll  motion 
and  one  for  pitch  motion)  and  a single  linear  spring  (for  heave 
motion,  i.e.,  vertical  translation  of  the  center  of  mass)  (Refs  3 
and  21).  A slightly  more  complex  model  has  also  been  used  which 
assumed  nonlinear  springs  (Ref  17:86).  The  most  elaborate  attempts 
to  model  the  dynamics  of  an  air  cushion  device  are  those  which 
mathematically  describe  the  pressure  distribution  and  air  flow  in 
and  around  the  trunk.  Reference  5 , pages  3 to  13,  presents  the 
salient  features  of  several  attempts  to  describe  the  trunk  pressure 
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and  air  flow  mathematically  in  order  to  predict  the  dynamics  of  the 
air  cushion  system.  Reference  5 also  presents  a sophisticated  ACR8 
model  which  takes  into  account  deformation  of  the  trunk  as  well  as 
air  supply  system  characteristics.  As  in  many  other  modeling 
problems,  the  more  complex  models  tend  to  have  enhanced  predictive 
capability  at  the  expense  of  time  and  effort  required  to  generate 
the  model. 

While  there  often  is  need  for  a precise  model,  there  is  good 
motivation  for  using  the  simplest  of  models,  the  linear  spring/ 
damper  model,  in  certain  applications.  The  slideout  of  an  ACRS- 
equipped  vehicle  is  claimed  to  be  one  such  application  for  the 
following  reasons.  First,  the  rate  of  change  of  vehicle  orientation 
is  typically  small  throughout  the  slideout,  except  during  and  just 
after  touchdown.  Certain  external  conditions  could  produce  a rapid 
rate  of  change  of  the  vehicle  orientation  (e.g. , severe  gusts  or 
irregular  landing  surfaces),  and  these  conditions  should  not  be 
analyzed  with  a linear  model.  Secondly,  if  a control  system  is  to 
be  designed  for  an  ACRS-equipped  vehicle,  the  most  dominant  charac- 
teristics of  the  ACRS  are  probably  sufficient  to  design  the  control 
system.  Once  the  control  system  is  designed,  the  control  system 
could  be  used  in  a simulation  with  a more  elaborate  trunk  model  in 
order  to  verify  the  design. 

Due  to  the  limited  amount  of  time  available  for  conducting 
research  for  this  report,  the  scope  of  the  research  has  been  limited 
to  slideouts  beginning  at  a point  just  after  touchdown  transients 
have  damped  out.  This  portion  of  the  slideout  allows  the  use  of  a 
linear  spring-damper  model  in  the  slideout  simulation. 
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The  remainder  of  this  chapter  consists  of  the  development  of  a 
linear  trunk  model  which  is  combined  with  the  aerodynamic  model  in 
Chapter  V in  order  to  analyze  the  performance  of  the  Jindivik  in 

I 

the  slideout  phase.  A description  of  some  experimentally  obtained 
data  is  the  starting  point  for  building  the  model,  since  the  dynamic 
behavior  of  the  model  should  resemble  that  of  the  actual  system. 

After  presenting  this  data  some  pertinent  geometry  is  introduced, 
followed  by  an  explanation  of  the  placement  of  the  spring-damper 
units  in  the  model.  A description  of  the  forces  produced  by  spring 
displacement,  displacement  rate,  and  sliding  friction  is  then  given. 

Finally,  the  model  dynamics  are  shown  to  be  essentially  the  same  as 
those  of  the  actual  ACRS-vehicle  combination  in  heave,  pitch,  and 
roll  modes. 

Experimentally  Obtained  Trunk  Data 

Unpublished,  experimentally  obtained  trunk  data  was  made 
available  for  this  report  by  the  Air  Cushion  Group,  AFFDL.  The 
methods  used  to  obtain  this  data  are  described  in  Reference  21. 

As  a value  for  the  heave  damping  coefficient,  Cy,  was  not  available, 
the  value  of  C^  for  a similar  trunk  was  used. 

In  addition  to  the  data  in  Table  III,  representative  values 
for  the  coefficient  of  friction,  u,  were  supplied  by  the  Air  Cushion 
Group/AFFDL.  A value  of  n = .3  was  recommended  for  the  forward  1/3 
of  the  trunk  since  this  portion  of  the  trunk  is  partially  supported 
by  an  air  cushion.  The  purpose  of  the  air  cushion  is  to  reduce 
friction  forward  of  the  center  of  mass.  By  shifting  the  frictional 
force  towards  the  aft  portion  of  the  trunk,  a more  directionally 


stable  situation  is  obtained  in  slideouts.  This  increase  in  stability 
is  derived  from  restoring  yaw  moments  which  are  produced  by  frictional 
forces  aft  of  the  center  of  mass  whenever  the  velocity  vector  of  the 
vehicle  is  not  aligned  with  the  longitudinal  axis  during  ground 
contact.  A value  of  u = .8  was  recommended  for  the  aft  2/3  portion 
of  the  trunk. 

Table  III 

EXPERIMENTALLY  OBTAINED  DATA  FOR  THE  JINDIVIX  ACRS  TRUNK  NUMBER  2 


% 

Heave  Stiffness  Coefficient 

12,300  lb/ft 

Ch 

Heave  Damping  Coefficient* 

480  lb-sec/ft 

Kp 

Pitch  Stiffness  Coefficient 

1252  ft-lb/deg 

CP 

Pitch  Damping  Coefficient 

24.7  ft-lb-sec/deg 

XR 

Roll  Stiffness  Coefficient 

90.4  ft-lb/deg 

CR 

Roll  Damping  Coefficient 

8.65  ft-lb-sec/deg 

*Value  was  not  calculated  for  ACRS  No.  2.  This  value  is  taken 
as  a representative  value  based  on  data  given  in  Ref  21. 

Description  of  the  Trunk 


The  trunk  is  constructed  of  neoprene  coated  nylon  fabric  with 
tire  tread  rubber  3/8  inch  thick  along  the  entire  underside  (Ref  17). 
Nozzles  are  drilled  through  the  brake  tread  in  the  forward  1/3  of 
the  trunk  to  provide  the  air  cushion  discussed  above.  The  trunk  is 
mounted  on  the  Jindivik  such  that  the  equilibrium  pitch  angle  is  0°. 

The  pertinent  dimensions  of  the  trunk  given  in  Figure  32  are 
obtained  from  Reference  18.  The  curve  formed  by  the  center  of  the 
brake  tread  is  used  throughout  the  remainder  of  this  report  for 


the  placement  of  all  spring-damper  units,  with  one  exception.  The 
exception  is  a spring-damper  'unit  that  supplies  a portion  of  the 
total  heave  force,  and  which  is  located  at  the  center  of  pressure. 


fuselage 


Center 

of 

Pressure 


NOTES 


1.  All  dimensions  in  feet 

2.  Sketch  on  right  not  to 
scale 


Fig.  32  Ge-v.: 


of  the  Jindivik  ACES  Number  2 


Placement  of  the  Spring-Damper  Units 

Although  there  are  many  ways  to  place  the  spring-damper  units 
in  order  to  represent  the  ACRS-vehicle  combination,  any  specific 
arrangement  should  have  the  following  particular  characteristics. 
First  and  foremost,  the  arrangement  must  exhibit  the  same  static 
and  dynamic  characteristics  as  the  actual  ACR3.  Secondly,  there 
should  be  a sufficient  number  of  units  in  order  t'  .'..present  ade- 
quately the  shape  of  the  brake-tread  center  line , ir.ce  distances 
from  the  center  of  mass  affect  the  torques  acting  on  the  vehicle. 
Finally,  a representative  number  of  spring-damper  nits  must  be 
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used  for  the  forward  1/3  of  the  tread  area  and  the  aft  2/3  of  tread 
area  in  order  to  account  for  differences  in  the  friction  coeffi- 
cients . 

Following  the  above  guidelines,  a total  of  14  spring-damper 
units  was  chosen.  In  order  to  meet  both  roll  and  pitch  stiffness 
and  damping  coefficients,  individual  units  around  the  periphery  of 
the  trunk  are  allowed  to  have  one  of  two  stiffness  and  damping 
coefficients.  Seven  of  the  units  are  located  near  the  extreme 
forward  and  extreme  aft  portions  of  the  trunk — these  spring-damper 
units  are  referred  to  as  "pitch"  springs  and  dampers.  Another  six 
spring-damper  units  are  located  primarily  along  the  sides  of  the 
trunk  and  are  referred  to  as  'holl'  units . The  fourteenth  unit  is 
located  at  the  static  center  of  pressure  of  the  trunk  (that  point 
through  which  all  of  the  force  produced  by  pressure  is  assumed  to 
act).  The  stiffness  and  damping  of  this  fourteenth  unit,  referred 
to  as  the  heave  spring-damper  unit,  are  chosen  so  as  to  bring  the 
total  stiffness  and  damping  of  all  units  to  the  correct  values  for 
the  heave  mode.  Figure  33  describes  a placement  which  satisfies 
all  the  requirements  stated  above,  and  this  arrangement  is  used  for 
all  slideout  simulations  in  this  report.  The  labeling  of  the  spring 
damper  units  is  as  follows:  1)  peripheral  spring-darrper  units  are 
given  increasing  numbers  from  fore  to  aft,  2)  a P or  an  S is  added 
to  the  number  to  identify  port  or  starboard  units,  3)  the  heave  unit 
is  designated  as  unit  number  88.  The  exact  locations  of  these 
spring-damper  units  are  given  in  Table  IV. 

Values  for  stiffness  coefficients  and  damping  coefficients 
are  obtained  in  the  following  manner.  If  angular  deflections  are 
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assumed  to  be  small,  the  following  set  of  equations  describes  the 
relationships  between  stiffness  coefficients  of  the  pitch,  roll,  and 


heave  springs  (kp,  kr,  and  kh,  respectively): 

2 


2 

kp  Exi 


kD  Ey^_ 


pitch  +kr£*i 
springs 

pitch 

springs 


roll  + khX§8  = Kp 
springs 


roll 

springs 


= Kr 


(77) 


(78) 


where  Kp  and  Kr  are  given  in  Table  III. 

In  a similar  manner,  the  relationships  between  pitch  damping 
coefficient,  roll  damping  coefficient,  and  the  individual  spring 


damping  coefficients  can  be  written 
2 


cD  Exi 


cD  zyi 


pitch  + Cr  EXi 
springs 

pitch  + °r  lyi 
springs 


roll  + cp?C88  CP 
springs 


roll 

springs 


= Cr 


(79) 

(80) 


where  Cp,  cr,  and  cr  represent  damping  coefficients  for  the  pitch, 
roll,  and  heave  dampers,  respectively.  Two  additional  equations 
express  the  relationships  between  the  stiffness  coefficients  and 
damping  coefficients  of  the  springs: 

7kp  + 6kr  + kh  = Kh  (81) 

7cp  + 6cr  + Ch  = Cr  (82) 

An  additional  constraint  arises  from  the  requirement  that  the  center 
of  pressure  (the  point  through  which  all  reaction  force  Is  assumed 
to  act)  lie  0.5  ft  forward  of  the  center  of  mass.  Mathematically 
this  requirement  is  expressed  as 


xcdkH  = hp  Exi 


pitch  + kr  CX1 
springs 


roll  + khX88 
springs 


(83) 
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Considering  equations  77  through  83  it  is  noted  that  there  are 
seven  linear  algebraic  equations  with  six  unknowns:  kp,  kr,  %, 

Cp,  cr,  and  c^.  Equations  77  through  82  lead  to  a unique  solution 
for  these  variables,  but  the  solution  does  not  necessarily  satisfy 
equation  83.  Therefore,  an  additional  variable  is  added.  The 
variable  chosen  is  the  x-distance  for  one  pair  of  roll  spring-damper 
units,  the  4P  and  4s  units.  The  values  for  these  two  spring-damper 
■units  in  Table  IV  are  a result  of  solving  the  seven  equations  (77 
through  83).  Solutions  for  the  six  variables  kp,  kr,  kft,  cp,  cr, 
and  ch.  are  given  in  Table  V. 

Table  V 

VALUES  OF  STIFFNESS  AND  DAMPING  COEFFICIENTS  FOR  PITCH, 

ROLL,  AND  HEAVE  SPRING-DAMPER  UNITS 


kp 

463.53  Ib/ft 

cp 

2.96  “>^ec 

kr 

264.75  lb/ft 

cr 

w,  lb-sec 
ft 

^h 

7466.8  lb/ft 

Ch 

264.5  lb^ec 

Reaction  and  Frictional  Forces  Acting  on  the  Trunk 

In  developing  the  reaction  and  frictional  forces  on  the  trunk, 
two  assumptions  are  necessary  concerning  the  spring-damper  units. 
First,  the  units  are  assumed  to  be  rigid  in  all  directions  except 
the  direction  perpendicular  to  the  surface  of  the  earth.  This 
assumption  makes  sense  physically  since  the  flexible  structure  of 
the  ACRS  is  such  that  the  pressure  distribution  under  a portion  of 
the  brake  tread  acts  normal  to  the  landing  surface.  A second 
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assumption  is  that  no  force  is  created  by  a particular  spring-damper 
unit  whenever  the  spring  displacement  alone  results  in  a negative 
force.  This  situation  is  interpreted  to  mean  that  the  spring  is  no 
longer  touching  the  ground  and,  hence,  is  incapable  of  transmitting 
any  force  to  the  vehicle. 

The  method  used  to  obtain  trunk  forces  and  moments  acting  on 
the  vehicle  consists  of  three  steps:  1)  transformation  of  body 
position,  velocity,  and  angular  rates  into  spring  displacements, 
displacement  rates,  and  Xg  - ye  velocities,  2)  calculation  of  spring 
reaction  and  frictional  forces  from  the  information  obtained  in  the 
first  step,  and  3)  transformation  of  the  forces  and  moments  acting 
on  the  vehicle  into  the  body  coordinate  system  from  the  information 
obtained  in  the  second  step.  Mathematical  relationships  needed  for 
each  of  these  steps  are  developed  below  for  a typical  spring-damper 
unit. 

Step  1.  Let  Ro  represent 
the  position  vector  from  some 
arbitrary'  point  on  the  surface 
of  the  (assumed  flat)  earth  to 
the  center  of  mass  of  the  vehi- 
cle. Further,  let  represent 
the  position  vector  from  the 
center  of  mass  of  the  vehicle 
to  the  attachment  point  of  any 
one  spring-damper  unit.  Tne 
sum  of  these  two  position  vec- 
tors is  itself  a position 


Fig.  3^  Position  Vector  from  an 
Arbitrary  Point  on  the  Surface  of 
the  Earth  to  a Spring-Damper  Unit 
Point  of  Attachment 
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vector,  r±,  from  the  arbitrary  point  in  the  Xg  - ye  plane  to  the 
point  on  the  body.  If  this  vector  is  expressed  in  the  earth  frame, 
the  third  (or  ze)  component  is  the  coordinate  of  the  attachment 
point  of  the  spring.  Figure  34  is  a sketch  representing  this 
situation,  and  the  following  mathematical  relationships  summarize 


the  preceding  comments: 

r±  = R0  + rb 


(84) 


where 


and 


Letting 


then 


rx 

*e 

xb 

ry 

= 

ye 

y b 

rz 

ze 

zb 

(85) 


The  only  component  of  interest,  however,  is  the  rz  component,  since 
this  component  is  necessary  to  determine  the  spring  deflection. 
Hence,  rz  is  expressed  as 


rz  = zg  + c31xb  + c32yb  + c33zb  (86) 

where  c31,  c32,  and  c33  are  the  third  row  elements  of  the  trans- 
formation matrix,  C^.  All  terms  except  Zb  on  the  RHS  of  equation 
86  are  known  quantities  since  ze  and  are  determined  during  each 
iteration  of  a simulation  run,  and  xb  and  yb  are  given  in  Table  TV. 
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A value  for  zb  is  arbitrarily  chosen  as  zero  feet;  hence,  the  spring- 
damper  units  are  assumed  to  be  attached  to  the  vehicle  at  the  x^,  - yb 
plane . 


J 


The  displacement  rates  of  the  springs  are  also  required  and 
these  can  be  obtained  by  taking  the  time  derivative  of  equation  84: 

/"9  D _ irtU 

(87) 


Letting 


and 


dt  ~ dt 

^ + 

m- 

r.  " 

rx 

h 

m- 

q 

ye 

dt 


then  the  last  tern  in  equation  87  is  the  only  term  left  to  be 
evaluated.  Ihe  differentiation  of  rb  must  be  carried  out  in  the 
(assumed  inertial)  earth  frame,  and  this  is  most  expeditiously 
accomplished  by  use  of  the  theorem  of  Coriolis  (Ref  23:437,  438): 

where  W is  the  angular  velocity  vector  of  the  body  frame  with  respect 
to  the  earth  frame.  Ihe  first  term  on  the  RHS  of  88  is  zero  since 
the  position  vector  in  the  body  frame  is  constant.  Ihe  second  term 
is  now  determined  by  recalling  from  Chapter  II  that 


W|b  = 


Then  the  cross  product  becomes 

[w  x rbjb  = 


Qzb  - Ryb 

Rxb  - Pzb 

Pyb  - Q*b 


(89) 


93 


which,  when  pre-mltiplied  by  C®,  is  in  a form  which  allows  an 
evaluation  of  equation  87: 


rx  xe  p Qzb  - Ryb 

ry  = Ye  + Ch  ^b  “ Pzb 

f z ze  D Pyb  - Qxb 


The  displacement  rate  of  any  specific  spring-damper  unit  is  given  by 
equation  90  since  all  variables  on  the  RHS  of  equation  90  are 
available  in  each  iteration  of  the  computer  simulation. 

Step  2,  The  forces  produced  by  the  deflection  and  the  deflec- 
tion rate  of  a particular  spring  are  easily  calculated  as  follows: 

Fd  = -k(rz  - rZQ)  + ^ (91 ) 


Fr  ~ ~crz 


where 


Fd  - force  due  to  spring  deflection 
Fr  - force  due  to  damper  deflection  rate 
rZo  - coordinate  of  spring  attachment  point  with  vehicle  at 
rest  on  the  ground 

Fdo  ~ spring  force  with  vehicle  at  rest  on  the  ground 
The  spring-damper  forces  are  expressed  in  the  earth  frame. 

The  only  remaining  force  to  consider  is  the  frictional  force. 
The  total  reaction  force,  Fd  + Fr>  snd  the  direction  of  the  velocity 
of  the  spring-damper  point  of  contact  are  necessary  in  the  computa- 
tion of  this  force.  The  components  of  velocity  in  the  xe  - ye  plane 
are  given  in  equation  90.-  Hence,  the  direction  of  the  velocity  in 
this  plane  is  conveniently  determined  by 


Y = tan' 


rife 
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where  y is  the  angle  between  the  velocity  vector  and  the  xe  axis 
with  a rotation  from  the  xe  to  the  ye  axis  being  the  positive  direc- 
tion. The  components  of  the  frictional  force  on  a given  spring 


damper  unit  are  given  by 


Ff x = -u  (Fd  + Fr  )cosy 
Ffy  = -u(Fd  + Fr)siny 


where 


Ff  , Ff  - components  of  frictional  force  in  the  xe  and  ye 
■X  y 

directions , respectively 

y - coefficient  of  friction  associated  with  a specific 
spring-damper  unit 

Step  3.  The  reaction  and  frictional  forces  derived  in  step  2 
are  all  expressed  in  the  earth  reference  frame.  The  purpose  of  step  3 
is  eventually  to  transform  these  forces  into  forces  and  moments  in  the 
body  reference  frame. 

The  forces  in  step  2 are  first  expressed  as  earth-referenced 
coordinates  of  a vector,  Ff,  as  follows: 


w - 


Fr  + Fd 


Transformation  of  these  forces  into  the  body  frame  is  easily  accom- 


plished as  follows: 


Xf  Ffx 

b 5 Yf  - Ce  FW 


e xy 
Fr  + Fd 
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The  moments  about  the  center  of  mass  of  the  vehicle  are  com- 
puted by  forming  the  cross  product  of  the  position  vector  to  the 
point  where  each  spring  contacts  the  ground,  and  the  force  vector: 


M = r*t  x ff 


(98) 


where 

**t  = 1)  + ce 

The  term  r^3  is  the  third  component  of  the  vector  r±  in  equation  84 
and,  furthermore j this  value  in  earth-frame  coordinates  is  available 
during  each  computer  iteration.  Expansion  of  equation  98  yields 

fytZf  - ztYf") 


o 

o 

"ri5 


(99) 


ztxf  - xtZf 


(100) 


ytYf  - ytxf 


where  Xf,  Yf,  and  Zf  are  the  coordinates  of  Ff  in  the  body  frame  (see 
equation  97  )>  and  xt,  yt>  and  zt  are  the  coordinates  of  r*t  in  the 
body  frame. 

This  concludes  the  three  steps  involved  in  calculating  the 
trunk  reaction  and  frictional  forces  in  the  body  frame.  The  FORTRAN 
implementation  of  these  results  is  given  in  Appendix  B. 


Verification  of  Trunk  Cynamics 

The  three  individual  dynamic  modes  for  which  data  is  available 
are  heave,  pitch,  and  roll,  as  given  in  Table  III.  A verification 
of  the  computer  program  with  the  trunk  model  included  is  desired  to 
ensure  that  no  programming  errors  have  been  conmitted.  Computer 
simulations  of  each  dynamic  mode  are  presented  in  this  section, 
which  verify  that  the  trunk  model  accurately  reflects  the  given 
heave,  pitch,  and  roll  data.  The  data  is  first  used  to  predict  the 
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damping  coefficient,  4 , and  period,  T,  for  each  mode,  after  which  a 
comparison  is  made  between  these  predictions  and  the  values  obtained 
by  computer  simulation. 

Prediction  of  heave  mode  dynamics.  In  the  heave  mode,  the  only 
degree  of  freedom  is  vertical  motion  of  the  center  of  mass.  A 


differential  equation  which  describes  this  mode  is 

mh  = -KHh  - Cyh  (101) 

where  h represents  the  deflection  of  the  spring  attachment  point  from 
its  equilibrium  position  and  the  positive  direction  of  h is  the 
opposite  of  ze.  A rearrangement  of  equation  101  yields 


which  has  the  form 


where 


••  Ck  . Kh 

h + — — - h + — h = 0 


h + 2SHUnHh  + = 0 


(102) 


(103) 


- heave  damping  ratio 
ojnpj  ~ heave  natural  frequency 

Equation  102  is  readily  solved  for  % and  Tr  (the  period  of  the  heave 


mode  oscillation) : 


2/mKH 

Tr  = — = i 


(104) 


(105) 


<1  - ^H2) 


Inserting  the  values  for  each  of  the  parameters  in  equations  104  and 
105,  the  following  experimentally-derived  values  for  4r  and  Tr  are 


obtained: 


4h  = 0.24 
Th  = 0.53  sec 
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Prediction  of  pitch  mode  dynamics.  The  equation  for  the  pitch 
mode  dynamics  is  similar  to  equation  101: 

IyyG  = -Kp0  - CP5  (106) 

Equation  106  is  solved  in  the  same  manner  as  equation  101  to  yield 

£p  = .062 
Tp  = .995  sec 

Prediction  of  roll  mode  dynamics.  The  equation  for  the  roll 
mode  dynamics  is  similar  to  equations  101  and  106: 

Ixx*  = "Kr$  " cr^  (!°7) 

Solving  this  equation  for  the  damping  ratio  and  period,  the  following 
values  are  obtained: 

Sr  = -102 
Tr  = 2.98  sec 

Simulation  results.  Computer  simulation  results  of  the  model 
in  each  dynamic  mode  are  presented  in  Figures  35,  36,  and  37.  These 
plots  are  treated  as  underdamped,  simple,  second-order  system  re- 
sponses and  are  graphically  evaluated  for  the  period,  T'f  and  the 
damping  ratio,  s',  which  are  compared  with  the  predicted  values 
computed  above.  The  manner  in  which  this  evaluation  is  performed  is 
now  described. 

With  the  assumption  that  the  time  responses  in  Figures  35  to 
37  represent  the  outputs  of  unforced,  second-order  systems,  the 
responses  are  of  the  form 

C(t)  = Ae-^  sin(o)(jt  + 4>)  (108) 
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where 


C(t)  - the  system  output 

A - a constant  determined  by  system  parameters 

. (108a) 

wd  - the  damped  natural  frequency  (aid  = ton*  1 - 52) 

$ - the  phase  angle 

The  product  of  and  the  sine  function  results  in  peak  values 

whenever  the  value  of  the  sine  function  is  near  ±1  and  the  damping 
ratio  is  small.  If  any  two  positive  peak  values  are  taken  from  the 
plots,  say  CCtj)  = Cj  and  C(t2)  = C2,  then  the  following  approximate 
equation  holds: 

£2.  Ae-Co>n(t2  - tj)  (109) 

An  expression  for  the  damping  ratio  is  determined  by  using  this 
relationship  and  the  expression  given  above  for 


C = 


ln(Cl/C2) 


V 


ln(Ci/C2)' 


4ir2(t2  - tl)2 
T 


(no) 


This  expression  is  used  to  evaluate  c sevei_J.  times  for  each  simu- 
lation, after  which  an  average  value,  t',  is  computed  in  order  to 
compare  the  predicted  value  and  the  simulation  value.  Furthermore, 
several  values  of  the  period  computed  for  each  simulation  were 
averaged  to  obtain  an  average  value,  T',  which  is  compared  with  the 
predicted  period,  T.  Thble  VI  lists  the  predicted  versus  the  simu- 
lation values  of  period  and  danping  ratio  for  each  dynamic  mode. 

Ihe  agreement  is  sufficiently  close  (less  than  5%  deviation  in  each 
case),  to  lead  to  the  conclusion  that  there  are  no  major  programming 
errors.  Therefore,  the  model  trunk  reaction  forces  are  essentially 
the  same  as  those  of  the  actual  ACRS- vehicle  combination. 
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Table  VI 


COMPARISON  OF  PREDICTED  AND  SIMULATED  DAMPING  RATIOS  AND  PERIODS 


Dynamic 

Mode 

Predicted 

Damping 

Ratio 

Simulation 

Damping 

Ratio 

Predicted 

Period 

(sec) 

Simulation 

Period 

(sec) 

HEAVE 

.2k 

.2k 

.53 

.53 

PITCH 

.062 

.060 

.995 

1.03 

ROLL 

.102 

.102 

2.98 

2.98 

Conclusion 

A reasonable  model  of  the  trunk  reaction  and  frictional  forces 
has  been  developed  which  is  claimed  to  be  realistic  for  small  angles. 
Since  test  results  in  Table  VT  are  based  on  peak  values  of  roll  and 
pitch  of  approximately  5°  and  3°,  respectively,  simulation  results 
are  not  guaranteed  to  be  accurate  outside  of  these  ranges. 

Analysis  of  the  dynamic  characteristics  of  the  Jindivik  in  the 
slideout  phase  is  now  possible.  Chapter  V uses  the  aerodynamic 
model  in  conjunction  with  the  trunk  model  to  perform  such  an  analysis. 
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V.  ANALYSIS  OF  THE  SLIDEOUT 


The  aerodynamic  model  developed  in  Chapter  II  and 
the  trunk  model  developed  in  Chapter  IV  are  combined  in 
this  chapter  for  the  purpose  of  analyzing  the  slideout 
stability  of  an  ACRS-equipped  aircraft.  Simulations  of 
the  Jindivik  aircraft  during  a slideout  are  analyzed 
which  serve  to  illustrate  a basic  stability  problem. 

Some  general  comments  concerning  the  nature  of  this 
problem  are  made,  followed  by  a proposed  solution. 

Finally,  simulation  results  are  presented  which  serve 
to  verify  the  effectiveness  of  the  proposed  solution. 

Analysis  of  the  ACF-S-Equlpped  Vehicle  in  the  Slideout 

Two  simulations  were  initially  performed  in  order  to  identify 
instabilities  in  the  slideout  phase.  Simulation  results  presented  in 
this  section  we re  obtained  with  the  configuration  and  initial  condi- 
tions given  in  Table  VII.  The  pitch  and  roll  control  systems  were 
omitted  from  the  first  simulation  in  order  to  establish  a baseline 
performance  capability  of  the  landing  system.  The  existing  pitch  and 
roll  control  systems  were  included  in  the  second  simulation  to  observe 
the  adequacy  of  these  controllers  in  the  slideout. 

Bare  Airframe  Slideout  Simulation 

Results  of  the  bare  airframe  simulation  are  presented  in  Figures 
38  through  U2.  The  overall  motion  of  the  vehicle  is  best  described 


Table  VII 

CONFIGURATION  AND  INITIAL  CONDITIONS  FOR  SLIDEOUT  SIMULATIONS 


Forward  Velocity 

200  ft/sec 

CG  Height1 

2.4977  ft 

Elevator  Deflection  Angle1 

-1.1354° 

Angle  of  Attack1 

-2.1302° 

Pitch  Angle1 

-2.1302° 

Thrust 

300.  lb 

Flap  Setting 

UP 

Weight 

2620  lb 

Friction  Coefficient 

(Forward  )2 

0.3 

Friction  Coefficient 

(Aft)2 

0.8 

lrThese  values  were  computed  using;  the  "STEADY  STATE"  feature 
of  the  computer  program,  EASY  (see  Appendix  B). 

2These  values  were  increased  for  the  individual  springs 
since  a significant  portion  of  the  total  reaction  force 
is  supported  by  the  spring  at  the  center  of  pressure, 
which  does  not  contribute  to  the  frictional  force. 
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by  separately  considering  the  longitudinal-related  dynamics  and  the 
lateral-directional  related  dynamics. 

Longitudinal  Dynamics.  Figure  38  depicts  the  time  histories  of 
the  xb  velocity,  U,  and  angle  of  attack,  a.  The  xt,  velocity  is  seen 
to  decrease  at  a fairly  constant  rate  from  the  initial  value  of 
200  ft/sec  to  a value  of  50  ft/sec,  where  the  simulation  is  ter- 
minated. The  average  deceleration  is  approximately  16.67  ft/sec 
(.52  g). 

The  angle  of  attack  gradually  decreases  during  the  first  6.H 
seconds  for  the  following  reason.  As  the  velocity  decreases,  the 
pitch-up  aerodynamic  moment  decreases,  while  the  nose-down  moment 
produced  by  sliding  frictional  forces  remains  relatively  constant. 

After  6.5  seconds,  large  fluctuations  in  a occur  due  to  roll  angle 
variations  which  are  discussed  in  the  section  on  lateral-directional 
dynamics.  Neither  the  deceleration  of  the  vehicle  nor  the  angle  of 
attack  variation  with  time  are  objectionable,  except  for  the  divergent 
tendency  of  the  angle  of  attack  after  five  seconds. 

The  pitch  angle  and  center  of  gravity  (eg)  height  are  presented 
in  Figure  39-  A steady  decrease  in  pitch  angle  results  for  the  same 
reason  that  the  angle  of  attack  decreases.  The  final  decrease  in  pitch 
angle  occurs  due  to  a large  deviation  in  roll  angle  which  is  discussed 
in  the  next  section. 

Although  the  pitch  angle  is  not  less  than  -3°  until  the  roll 
angle  becomes  quite  large,  the  nose-down  attitude  is  cause  for  concern, 
since  only  a slight  perturbation  could  result  in  the  nose  of  the  air- 
craft contacting  the  runway.  A control  system  could  be  designed  to 
maintain  a desired  pitch  attitude  during  the  initial  portion  of  the 
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Fig.  38  Forward  Velocity  and  Angle  of  Attack  Time  Histories 
During  Slideout  with  no  Control  System 


Fig.  39  Pitch  Angle  and  Height  of  CG  Time  Histories  During 
Slideout  with  No  Control  System 


slideout,  but  the  largest  negative  pitch  deviations  occur  toward  the 
end  of  the  slideout  where  there  is  insufficient  dynamic  pressure  for 
the  elevator  to  be  effective.  Hence,  no  specific  design  is  developed 
for  controlling  the  pitch  angle  in  this  report. 

The  eg  height , h^,  is  initially  lower  than  the  static  eauili- 
brium  height  of  2.51  ft  due  to  aerodynamic  forces.  The  springs  of 
the  trunk  model  are  initially  compressed,  but  they  extend  toward 
their  equilibrium  values  as  the  aerodynamic  forces  lessen.  The  final 
large  increase  in  eg  height  is  a result  of  the  large  roll  angle  dis- 
cussed in  the  next  section. 

Lateral-Directional  Dynamics 

A small  (-3  ft/sec)  crosswind  gust  was  included  in  the  simula- 
tions in  order  to  excite  the  unstable  characteristics  of  the  ACRS- 
equipped  vehicle.  The  gust  was  applied  two  seconds  into  the  simula- 
tion when  the  forward  velocity  was  approximately  160  ft/sec.  Figure 
40  depicts  time  histories  of  the  roll  angle,  $,  and  yaw  angle,  V. 

The  roll  angle  is  seen  to  remain  at  0°  until  the  gust  Is  applied. 
Since  the  gust  produces  an  instantaneous  increase  in  the  sideslip 
angle,  a negative  rolling  moment  is  also  produced  which  causes  the 
aircraft  to  roll  to  the  left.  The  aerodynamic  characteristics. of  the 
Jindivik  are  such  that  an  increase  in  3 produces  a negative  yawing 
moment.  The  negative  roll  angle  increases  the  reaction  and  frictional 
forces  on  the  left  (pert)  side  of  the  vehicle  which  also  cause  a 
negative  yaw  moment  about  the  eg.  These  yaw  moments  produce  the  yaw 
rate,  R,  shown  in  Figure  4l,  which  results  in  the  negative  increase 
in  yaw  angle  in  Figure  40.  As  the  velocity  decreases  the  aerodynamic 
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yawing  moment  decreases,  while  the  yawing  effect  of  unbalanced 
frictional  forces  remains  relatively  constant.  Thus,  the  yaw  rate 
tends  to  follow  the  roll  angle  throughout  the  sliaeout  and  the  shape 
of  the  curves  for  $ and  R are  seen  to  be  very  similar. 

The  roll  angle  influences  yaw  rate,  but  the  roll  angle  itself 
is  influenced  by  both  lateral  frictional  forces  and  spring  reaction 
forces.  Thus,  as  the  yaw  angle  increases  negatively,  the  frictional 
forces  are  acting  to  produce  a positive  rolling  moment.  This  moment, 
along  with  trunk  reaction  forces  causes  the  roll  angle  to  increase 
from  its  negative  value  of  -1.5°  at  4 seconds,  to  7°  at  6.7  seconds. 
The  reaction  torque  becomes  larger  than  the  roll  torque  produced  by 
lateral  frictional  forces  just  after  6 seconds  and  the  roll  angle 
begins  to  decrease  sharply  at  this  point. 

Figure  42  depicts  the  lateral  displacement  of  the  vehicle,  ye, 
versus  distance  "down  the  runway",  xe.  It  is  seen  that  a sidewind 
gust  of  only  -3  ft/sec  produces  a deviation  of  -13  feet  from  the 
runway  centerline  (the  xe  axis). 

Airframe  with  Existing  Autopilot  Simulation 

Results  of  the  simulation  of  the  airframe  and  existing  autopilot 
are  presented  in  Figures  43  through  47.  The  longitudinal  dynamics 
are  very  similar  to  those  of  the  bare  airframe,  while  there  are  signi- 
ficant differences  in  the  lateral-directional  dynamics. 

Longitudinal  Dynamics.  Figure  43  presents  the  same  variables  as 
Figure  38  — U and  a . These  two  figures  may  be  compared  to  see  that 
the  trend  is  the  same  for  both  variables.  The  few  differences  which 
do  exist  are  primarily  a result  of  smaller  fluctuations  in  the  roll 
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Fig.  42  Lateral  Versus  Longitudinal  Displacement  During  Slideout 

with  No  Control  System 
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lgle  of  Attack  Time  Histories 
Existing  Autopilot 


Fig.  4*)  Pitch  Angle  and  Height  of  CG  Time  Histories  During 
Slideout  with  Existing  Autopilot 


angle.  These  differences  are,  1)  a more  linear  rate  of  change  In  U 
near  the  end  of  the  simulation,  and  2)  a more  stable  time  history  of 
a,  although  a is  still  somewhat  erratic  towards  the  end  of  the  simu- 
lation. 

The  pitch  angle  and  eg  height  depicted  in  Figure  44  are  seen  to 
fluctuate  considerably  less  with  the  pitch  and  roll  control  systems 
engaged  than  in  the  bare  airframe  simulation.  As  in  the  case  of  U 
and  a,  the  fluctuations  are  less  severe  primarily  because  the  roll 
angle  deviations  are  smaller. 

Lateral-Directional  Dynamics.  While  the  shape  of  the  roll  angle 
curve  in  Figure  45  is  similar  to  the  roll  angle  curve  in  Figure  40, 
the  magnitudes  are  considerably  smaller  in  Figure  45  due  to  the  action 
of  the  roll  control  system.  A dramatic  decrease  in  yaw  angle,  yaw 
rate,  and  sideslip  angle  deviations  is  also  noted  in  Figures  45  and  46. 

Figure  47  presents  the  lateral  versus  longitudinal  displacements 
of  the  vehicle  for  the  autopilot-engaged  simulation.  While  most 
other  variables  in  this  simulation  deviate  less  from  their  nominal 
values  than  they  do  in  the  bare  airframe  simulation,  the  final  lateral 
displacement  is  only  mildly  reduced  (from  -13  ft  to  -12  ft).  The 
primary  reason  for  the  small  difference  in  the  two  cases  is  that  the 
large  positive  yaw  angle  near  the  end  of  the  autopilot-off  simulation 
curtails  the  lateral  displacement  rate.  Due  to  this  large  roll  angle, 
however,  the  accuracy  of  the  simulation  is  questionable  in  this  area. 

Conclusions.  Two  primary  conclusions  are  drawn  from  these  simu- 
lations. First,  both  configurations  are  unacceptable  due  to  the 
divergent  tendency  of  the  roll  angle.  Although  in  the  latter  simu- 
lation the  roll  angle  achieves  a final  value  of  only  -6° , it  is 
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Fig.  46  Sideslip  Angle  and  Yaw  Rate  Time  Histories 
During  Slideout  with  Existing  Autopilot 


Fig.  l[7  lateral  Versus  Longitudinal  Displacement  During  Si ideout 

with  Existing  Autopilot 


decreasing  at  the  rate  of  -1.5  /sec.  Since  the  magnitude  of  the 
crosswind  gust  is  only  -3  ft/sec,  it  is  highly  probable  that  a larger 
gust  will  cause  the  roll  angle  to  become  unacceptably  large. 

A second  conclusion  is  that  the  roll  angle,  roll  rate,  yaw  angle, 
yaw  rate,  and  side  velocity  are  all  important  in  the  lateral-direc- 
tional dynamics.  The  descriptions  given  above  of  the  simulation 
results  necessarily  involve  each  of  these  variables  to  explain  the 
behavior  of  the  system.  Hence,  in  order  to  control  the  lateral- 
directional  dynamics  of  an  ACRS-equipped  vehicle,  these  variables 
must  all  be  taken  into  account. 

Analysis  of  the  Lateral-Directional  Instability 

In  order  to  gain  an  understanding  of  the  cause  of  the  instability 
discussed  in  the  previous  section,  a linearized  analysis  of  the 
slideout  was  performed.  The  analysis  was  carried  out  on  a digital 
computer  by  using  the  LINEAR  ANALYSIS  feature  of  EASY  (see  Appendix 
3) . The  computer  program  produced  a linearized  approximation  of  the 
state  model,  as  well  as  the  eigenvalues  of  this  linear  model,  at 
forward  velocities  of  200,  160,  120,  80,  and  40  ft/sec.  A sketch  of 
the  location  of  these  eigenvalues  in  the  complex  plane  is  presented 
in  Figure  48. 

From  this  sketch  the  following  observations  are  made:  1)  the 
eigenvalues  do  not  vary  greatly  with  velocity,  2)  three  eigenvalues 
are  near  the  origin,  and,  3)  one  complex  pair  of  eigenvalues  has  a 
positive  real  part.  Elaboration  of  each  of  these  observations  is  in 
order . 

The  first  observation  is  used  in  the  next  section  is  designing 
an  optimal  controller  for  the  landing  phase.  Since  the  variation  of 
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eigenvalues  is  small,  a single  set  of  constant  feedback  gains  can  be 
used  for  the  entire  slideout,  thus  greatly  reducing  the  complexity  of 
the  controller. 


The  three  eigenvalues  near  the  origin  are  associated  with  the 
variables  xe,  ye,  and  V.  This  conclusion  is  based  on  the  fact  that 
the  dynamics  of  the  system  are  virtually  independent  of  the  values  of 
these  variables.  Another  way  of  expressing  this  is  as  follows:  since 
Xe  and  ye  are  the  integrals  of  xe  and  ye,  respectively,  and  since  ¥ is 
approximately  the  integral  of  R,  a free  s would  appear  in  the  denomina- 
tion of  any  transfer  function  involving  xe,  ye,  or  ^ as  output  vari- 
ables. This  association  of  particular  variables  with  specific 
eigenvalues  is  used  in  the  next  section,  where  it  is  desired  to 
construct  a reduced  state  model  for  the  purpose  of  designing  a 
lateral-directional  control  system. 

The  third  observation  stated  above  is  the  most  important  result 
of  the  linearized  analysis  of  the  slideout.  The  pair  of  eigenvalues 
in  the  right-half  complex  plane  are  associated  with  the  instability 
noted  in  the  slideout  simulations.  The  location  of  these  eigenvalues 
must  be  modified  such  that  either  1)  there  are  no  eigenvalues  with 
positive  real  parts,  or  2)  the  magnitude  of  the  real  parts  of  any 
eigenvalues  in  the  right-half  complex  plane  are  sufficiently  small 
such  that  the  slideout  is  terminated  before  the  system  deviates 
excessively  from  the  desired  state.  The  optimal  control  design 
^ developed  later  in  this  chapter  effectively  relocated  these  eigenvalues 

into  the  left  half  of  the  complex  plane. 
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Design  of  a Lateral-Directional  Control  System  for  the  Slideout 


r 


Due  to  time  constraints  in  the  research  phase  of  this  report, 
only  one  design  technique  has  been  attempted.  While  the  magnitude 
of  the  problem  has  been  significantly  reduced  by  the  development  of 
a slideout  model,  much  work  remains  to  be  done  in  order  to  develop  a 
practical  controller.  The  optimal  control  design  presented  below  is 
only  one  approach  to  the  problem,  and  while  it  appears  to  be  highly 
successful,  other  design  approaches  as  well  as  variations  to  this 
design  should  be  pursued. 

Continuous  optimal  control  theory  is  used  in  designing  a lateral- 
directional  control  system.  As  this  theory  is  well  developed  in  the 
literature  (Refs  13  and  2 ),  only  a brief  description  of  this  method 
is  discussed. 

Problem  Formulation.  The  overall  objectives  to  be  met  in  the 
design  process  are  1)  to  keep  the  vehicle  within  the  confines  of  some 
(as  yet)  unspecified  landing  area,  and  2)  to  maintain  an  acceptable 
orientation  of  the  vehicle  during  the  entire  slideout.  It  is  assumed 
at  the  outset  that  a human  controller  on  the  ground  is  to  have  some 
control  of  the  vehicle.  Hence,  the  design  of  a control  scheme  need 
meet  only  the  second  objective  stated  above,  provided  a sufficient 
means  for  controlling  the  system  on  the  ground  is  available.  There- 
fore, direct  control  of  lateral  displacement,  ye,  is  not  attempted  in 
the  design  of  a control  system,  since  it  is  assumed  that  lateral 
displacement  can  be  controlled  by  controlling  the  roll  angle.  This 
appears  to  be  a reasonable  assumption  in  light  of  the  slideout  results 
presented  earlier,  but  this  assumption  should  be  verified. 


While  the  entire  system  could  be  treated  as  a multiple-input/ 
multiple-output  system,  it  is  possible  to  decompose  the  system  into 
longitudinal  and  lateral-directional  reduced-state  models.  Each  of 
these  systems  can  be  treated  as  a single-input/multiple-output 
system,  and  this  approach  is  taken  here. 

Four  longitudinal  states  were  selected  on  physical  grounds : 9 , 
Q,  V/,  and  hcg  (altitude  of  the  center  of  gravity).  These  four  states 
two  of  the  three  states  associated  with  eigenvalues  equal  to  zero 
(xe  and  ye),  and  U were  all  eliminated  from  the  full  12-state  system 
in  an  attempt  to  identify  the  lateral-directional  dominant  states. 

The  remaining  states  are  $,  4',  V,  P,  and  R.  Various  attempts  to 
reduce  the  lateral-directional  states  further  were  fruitless.  These 
five  states  may  well  be  the  minimum  number  of  states  which  can  be 
used  to  design  an  effective  control  system.  Fortunately  each  of  thes 
five  states  is  capable  of  being  sensed  onboard  the  vehicle,  although 
sensor  accuracy  must  be  considered.  A linear  analysis  was  performed 
using  the  digital  computer  program,  EASY  (see  Appendix  B),  at  U = 120 
ft/sec  in  order  to  verify  that  the  5 states  chosen  are  the  dominant 
lateral-directional  states.  The  linearized  state  equation  at  this 
velocity  is 


0. 

0.  0. 

1. 

0. 

0 

0. 

0.  0. 

0. 

1. 

0 

X = 

.0367 

.3600  -.0770 

.0052 

-2.090 

x + 

0 

-3.050 

-3,580  -.0675 

-1.110 

.7150 

-82. 

3.^33 

-1.030  -.0500 

.2784 

-.1572 

6.6 

where  xr 

r = [*, 

V,  V,  P,  R],  and 

u = 6A. 

u (110) 


The  eigenvalues  of  this  state  model  are  given  below,  along  with 


the  nearest  corresoonding  full-state  eigenvalues  for  comparison. 

5 State  Lateral-  Corresponding 

Directional  Eigenvalues  Full-State  Eigenvalues 


* -.09015  -.008596 

* .7M7±j  1.505  .6460±j  1.3^1 

-1.372±j 1.822  -1.220+J 1.808 


Although  the  agreement  is  not  excellent,  the  locations  of  the  5-state 
system  roughly  approximate  the  dominant  lateral-directional  dynamics. 

The  system  in  equation  110  is  in  a form  which  allows  optimal 
control  theory  to  be  applied.  Since  the  Immediate  objective  is  to 
regulate  the  system  (i.e.,  maintain  the  system  states  near  the  origin 
of  the  state  space),  the  problem  is  properly  termed  a linear,  time- 
invariant,  optimal  regulator  problem.  A mathematical  description  is 
given  by  first  defining  a quadratic  performance  index,  J,  for  the 
problem  at  hand: 


J 


[x-Qx  + u^]dt 


(111) 


where  Q is  a weighting  matrix.  The  solution  to  the  optimal  linear 
regulator  problem  with  this  performance  index  is  in  the  form  of  a 
control  law  with  a constant  feedback  gain  matrix , F : 


u*  = -Fx 


(112) 


where  the  symbol  is  used  to  denote  the  optimal  control.  A feed- 
back gain  matrix  was  determined  by  using  OPTCON,  a digital  computer 
program  (Ref  22).  For  reasons  discussed  later,  the  weighting  matrix 
v/as  chosen  to  be 


q = 


‘0.1  0000 
01000 
00100 
00010 
00001 
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which  resulted  in  the  feedback  matrix 


F = 


-3.6156 

-1.7124 

-.2919 

-1.2103 

-2.1646 


The  use  of  these  feedback  gains  results  in  a new  set  of  eigenvalues 

for  the  5-state  system  as  shown  below. 

Eigenvalues  of  5-State  Eigenvalues  of  5-State 
System  without  Control  System  with  Optimal  Control 


-.09015 

.7477+31.505 

-1.372+J1.822 


-.166 

-1.52 

-82.24 

-1.19+J1.85 


The  elements  of  the  weighting  matrix  were  chosen  such  that  all 
state  deviations  are  penalized  equally  except  for  roll  angle.  The 
choice  of  the  ’weighting  matrix  is  arbitrary  so  long  as  Q is  positive 
semi-definite,  but  the  performance  of  the  system  is  significantly 
affected  by  the  elements  of  Q.  VJhile  the  particular  values  chosen 
for  the  Q matrix  in  this  control  design  provide  a stable  system,  a 
more  refined  design  process  would  involve  a tuning  of  Q to  provide 
even  better  dynamic  characteristics. 


Verification  of  the  Lateral-Direct ional  Control  System 

Three  simulation  runs  were  performed  with  varying  crosswind 
magnitudes  introduced  two  seconds  into  the  simulation.  The  existing 
pitch  control  system,  and  the  lateral-directional  control  system 
designed  above  were  used.  P.esults  of  these  simulations  are  shown  in 
Figures  49  through  57.  The  longitudinal  states  are  not  shorn  since 
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Fig.  49  Roll  Angle  and  Yaw  Angle  Tjme  Histories 
During  Slideout  with  Optinal  Controller  (-3  ft/sec  Gust) 


SLIOEOUT  S IMULRTION  KITH  OPTIMAL  CONTROLLER 
ca/cn/n  17.13.16.  sir»unioN  oi5Puw 
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31  Lateral  Versus  Longitudinal  Displacement  During  Slideout 
with  Optimal  Controller  (-3  ft/sec  Gust) 


Pig.  53  Sideslip  Angle  and  Yaw  Rate  Time  Histories 
During  Slideout  with  Optical  Controller  (-10  ft/sec  Gust) 


Fig.  54  Lateral  Versus  Longitudinal  Displacement  During  Slideout 
with  Optimal  Controller  (-10  ft/sec  C-ust) 
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their  time  responses  are  very  similar  to  those  in  the  case  of  the 
slideout  with  the  existing  autopilot. 

The  dramatic  improvement  in  stability  is  readily  apparent. 

The  lateral  deviation,  ye,  is  -4.7  ft  with  a -3  ft/sec  constant 
crosswind  gust  at  2 seconds,  which  is  less  than  505  of  the  deviation 
achieved  without  the  optimal  controller.  The  deviations  are  consi- 
derably higher  for  crosswind  gusts  of  -10  ft/sec  and  -25  ft/sec 
(-20  ft  and  -44  ft,  respectively),  but  it  is  believed  that  the 
stability  of  the  vehicle  is  such  that  a human  operator  with  control 
of  the  roll  angle  could  effectively  control  the  ground  track. 

Although  the  magnitude  of  the  aileron  deflection  angle  and 
angular  rate  are  much  less  than  their  maximum  values  (±10°,  53°/sec, 
respectively ) during  the  entire  slideout  with  crosswind  values  of 
-3  ft/sec  and  -10'  ft/s.ec,  the  deflection  angle  did  exceed  its  limit 
in  the  -25  ft/sec  gust  simulation  at  a forward  velocity  of  approxi- 
mately 60  ft/sec.  Thus,  there  is  a need  for  low  speed  roll  control 
since  the  ailerons  become  ineffective  at  low  speeds. 

Low  speed  roll  control  can  be  achieved  through  the  use  of  roll 
thrusters.  The  output  of  these  devices  is  similar  to  the  yaw 
thruster  output  in  that  the  forces  are  constant  in  either  direction. 
A design  of  a low  speed  roll  control  system  is  not  undertaken  here, 
but  a design  procedure  very  similar  to  that  given  in  Chapter  III 
would  likely  be  very  fruitful. 

Conclusions 

The  following  conclusions  are  made  based  on  research  efforts 
presented  in  this  chapter: 
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1.  The  ACRS-equipped  Jindivik  is  unstable  during  slideout 
with  the  existing  control  system  engaged  or  disengaged. 

2.  A single  complex  pair  of  eigenvalues  lies  in  the  right 


half  complex  plane  which  is  a result  of  the  lateral- 
directional  state  variable  relationships. 

3.  A reduced-state  optimal  control  design  is  effective  in 
stabilizing  the  full-state  system. 

Some  method  of  roll  control  at  low  speed  is  required. 
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VI.  SUMMARY  Aid  RECOMMEldATIONS 
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In  this,  the  sixth  and  aonaluding  chapter,  highlights 
of  the  preceding  five  chapters  are  summarized,  and  areas 
for  further  research  are  recommended.  In  the  summary 
section , some  comments  on  the  generality  of  the  findings 
of  this  report  are  given.  The  recommendations  for  further 
research  are  given  in  two  sections : l)  recommendations 
for  additional  in-flight  study,  and  2)  recommendations  for 
further  investigation  of  the  slideout. 

Su mmarv 

Following  the  introductory  chapter,  which  provides  background 
inforrration  on  air  cushion  recovery  systems,  a nonlinear,  six 
degree  of  freedom  aerodynamic  model  is  derived.  The  aerodynamic 
model  includes  the  existing  autopilot  control  laws.  Furthermore, 
the  model  allows  selection  of  thrust  and  flaps  setting. 

Chapter  III  provides  1)  an  analysis  of  the  bare  airframe, 

2)  an  analysis  of  the  existing  autopilot,  3)  an  aerodynamic  in- 
stability, 4)  a proposed  solution  to  the  stability  problem,  and 
5)  a verification  of  the  proposed  solution.  Although  the  analyses 
in  Chapter  III  are  based  on  the  Jindivik  aircraft  in  particular,  the 
yaw  instability  discovered  is  primarily  attributable  to  the  addition 
of  an  ACRS.  Hence,  the  problem  Is  of  a general  nature,  and  is 
likely  to  occur  again  in  the  development  of  air  cushion  recovery 
systems . 
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The  fourth  chapter  provides  a spring-damper  approach  to  modeling 
ground  reaction  and  frictional  forces  which  are  transmitted  to  the 
ACES  trunk  during  a slideout.  A specific  trunk  model  for  the 
Jindivik  aircraft  is  achieved,  but  the  approach  taken  can  be  applied 
to  other  ACPS-equipped  vehicles.  While  experimentally  obtained  data 
is  used  in  Chapter  IV  to  model  the  trunk,  new  ACRS  designs  may  be 
analytically  evaluated  by  the  use  of  models  which  calculate  the 
pressure  distribution  and  deformation  of  the  trunk.  The  results  of 
these  simulations  could  provide  sufficient  information  on  the  roll, 
pitch,  and  heave  characteristics  of  the  ACPS-vehiele  combination  to 
allow  a spring-damper  model  to  be  developed. 

Chapter  V presents  a cursory  analysis  of  the  Jindivik  aircraft, 
equipped  with  an  ACRS,  during  the  slideout.  A stability  problem  is 
identified  by  both  simulations  of  the  slideout  and  by  analysis  of 
linearized  models  at  different  forward  velocities.  It  is  not  known 
at  this  time  whether  the  cause  of  this  instability  is  peculiar  to 
the  Jindivik,  or  if  this  is  a more  general  trait  of  air  cushion 
recovery  systems. 

Recommended  Additional  Investigation  (In-Flight) 

Additional  aerodynamic  data  on  the  Jindivik  equipped  with  an  air 
cushion  recovery  system  is  highly  recommended.  An  interesting  thesis 
for  an  aerodynamics  student  might  well  be  an  analytical  study  of  the 
effects  of  an  ACRS  on  the  aerodynamic  characteristics  of  an  aircraft. 
Wind  tunnel  testing  is  also  needed,  especially  for  the  lateral- 
directional  aerodynamic  coefficients, 

* 
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The  yaw  controller  design  in  Chapter  III  requires  considerably 
more  attention  than  could  be  given  in  this  report.  The  following 
areas  are  specifically  recommended  for  future  investigation: 

1.  Development  of  a control  law  for  the  yaw  thruster  when  the 
aircraft  is  in  a turn. 

2.  Determination  of  the  effects  of  varying  airspeed  on  the 
performance  of  the  yaw  controller. 

3.  Implementation  problems  which  might  arise,  such  as,  1) 
sensor  accuracy  required  versus  present  capabilities,  and 
2)  minimum  sampling  time  required  for  digital  implementa- 
tion. 

Recommended  Additional  Investigation  (Slideout) 

Many  challenging  and  stimulating  control  problems  arise  from 
the  development  of  the  aerodynamic  and  trunk  models.  Some  of  the 
more  important  problems  which  require  additional  investigation  are 
given  below. 

1.  Eigenvalue  sensitivity  to  parameter  variation. 

While  considerable  effort  has  been  expended  in  modeling 
the  aerodynamic  and  trunk  forces,  some  of  the  parameters  needed 
could  not  be  calculated  with  great  precision.  Since  the  loca- 
tion of  the  eigenvalues  is  instrumental  in  designing  a controller, 
a knowledge  of  how  the  eigenvalues  vary  with  certain  parameters 


1H0 


is  needed. 


2.  Sensitivity  of  optimal  controller  performance  to  eigenvalue 
variations . 

The  optimal  control  technique  described  in  Chapter  V may 
be  very  sensitive  to  errors  in  the  eigenvalue  locations.  ! 

Specifically,  if  the  optimal  controller  acts  so  as  to  cancel 
the  eigenvalues  in  the  right-half  complex  plane,  inexact 
cancellation  might  occur  and,  thus,  the  system  might  still  be 
unstable . 1 

3.  Variations  in  the  optimal  control  design. 

Several  variations  in  the  optimal  control  problem  formula- 
tion could  be  analyzed.  Among  the  most  promising  of  these  are 

(i)  Allowing  the  control  matrix,  B,  to  vary  with  time. 

(Since  the  lateral-directional  control  (aileron 
deflection  angle)  varies  with  the  square  of  total 
velocity,  the  use  of  a constant  B matrix  in  Chapter  V 
is  rather  unrealistic . ) 

(ii)  Changing  the  limits  on  the  quadratic  performance 

integral  (equation  111)  to  a finite  time.  (The  feed- 
back gains  would  become  time  varying,  but  improved 
performance  might  be  obtained.) 

(iii)  Tuning  the  weighting  matrix,  Q. 

4.  Evaluation  of  other  means  of  control  during  the  slideout. 

The  yaw  thruster  is  capable  of  delivering  a small  force 
(approximately  10  lbs)  and  a yaw  moment  (approximately  100  lbs) 
with  the  engine  at  idle.  This  force  and  moment  might  be  used 
effectively  at  lower  speeds. 
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Also,  if  a thrust  reverser  were  to  be  utilized,  the  engine 
could  be  operated  at  a higher  thrust  level,  and  a significant 
yawing  moment  would  be  available  from  the  yaw  thruster. 

While  the  ailerons  have  been  shown  to  be  effective  in 
controlling  the  vehicle  over  a wide  range  of  velocity,  roll 
control  is  needed  at  low  speeds.  The  use  of  roll  thrust  devices 
is  a recommended  area  for  further  investigation. 

5.  Analysis  of  the  landing  transient. 

While  the  trunk  model  described  in  this  report  is  not 
suitable  for  simulating  the  landing  transient,  a nonlinear 
formulation  of  spring  damping  and  stiffness  coefficients  might 
be  all  that  are  necessary  to  produce  a reasonable  model.  For 
example,  the  spring  forces  could  be  modeled  as  being  propor- 
tional to  the  cube  of  the  spring  displacement  (in  compression) . 

6.  Recommendations  for  pitch  control  technio.ues. 

As  mentioned  in  Chapter  V,  the  sliding  friction  force  pro- 
duces a large  nose-down  pitching  moment  during  slideout.  A 
successful  means  of  controlling  this  pitch  is  needed. 

Conclusion 

Clearly  there  are  many  problems  remaining  to  be  solved  before 
air  cushion  landing  gear  can  be  considered  as  an  effective  replace- 
ment for  conventional  gear.  The  problems,  however,  do  not  in  any 
way  appear  to  be  insurmountable. 
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Appendix  A 


Pertinent  Aerodynamic  Data 
for  the 

Jindivik  Drone  Aircraft 

This  appendix  consists  of  two  main  sections.  The  first  section 
contains  tabular  listings  of  the  wind  tunnel  data  and  curve  fits  to 
this  data.  The  second  section  presents  certain  aerodynamic  data 
provided  by  the  manufacturer  of  the  Jindivik,  as  well  as  some 
derived  aerodynamic  parameters. 
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All  wind  tunnel  data  was  taken  by  the  Government  Aircraft 

Factories  of  Australia  (Ref  15).  A one-quarter  scale  model  was  used. 

i 


The  vehicle  was  configured  with  an  ACRS  trunk  and  a 3-inch  (full 
scale)  pipe  on  the  port  side.  The  pipe  was  added  to  the  vehicle  to 
supply  bleed  air  from  the  turbine  engine  to  the  trunk  of  the  ACRS. 

Pertinent  wind  tunnel  data  is  presented  in  the  following 
tables.  Table  I-A  is  a listing  of  the  longitudinal  aerodynamic 
coefficients  (Cl,  Cq,  and  CM)  versus  angle  of  attack,  a,  for  a 
flaps-up  configuration.  Table  II- A is  a listing  of  the  same 
coefficients  but  with  the  flaps  down.  Variations  in  these  coeffi- 
cients with  sideslip  angle  are  available,  but  these  variations  are 
small  and  are  therefore  neglected. 

Tables  III-A  through  VI-A  consist  of  lateral-directional  aero- 
dynamic coefficients  versus  sideslip  angle,  8,  for  the  flaps-up 
configuration.  The  same  information  is  given  in  Tables  VII-A 
through  X-A  but  for  a flaps-down  configuration.  These  tables  show 
a variation  in  the  lateral-directional  coefficients  with  angle  of 
attack  for  both  flap  settings. 

Tables  XI-A  through  XIII-A  consist  of  changes  in  the  longi- 
tudinal aerodynamic  coefficients  with  elevator  deflection  angle,  6g, 
for  the  flaps-up  configuration.  The  flaps-down  configuration  is  not 
significantly  different  from  the  flaps-up  configuration  and,  hence, 
is  not  shown  here. 
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Table  I- A 

LONGITUDINAL  AERODYNAMIC  COEFFICIENTS  (Flaps-Up) 


ANGLE  OF 


:K  (DEGREES) 

CL 

CD 

CM 

• 09 

.1430 

.0940 

.0562 

2.19 

• 3143 

.1027 

.0415 

4.30 

.4982 

.1114 

.0367 

6.41 

.6770 

.1585 

■ 0326 

Table  II-A 

LONGITUDINAL  AERODYNAMIC  COEFFICIENTS  (Flaps-Down) 


ANGLE  OF 


OK  (DEGREES) 

cL 

CD 

Cm 

-1.75 

.4096 

.1609 

-.0063 

' -33 

.5443 

.1650 

-.0071 

2.45 

.7310 

.1884 

-.0122 

4.55 

.9107 

.2218 

-.0222 
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Table  xi-A 


CHANGE  IN  LIFT  COEFFICIENT  WITH  ELEVATOR  DEFLECTION  ANGLE 


Table  XII-A 

CHANGE  IN  DRAG  COEFFICIENT  WITH  ELEVATOR  DEFLECTION  ANGLE 


CHANGE  IN  DRAG  COEFFICIENT 


-.0183 

-.0201 

-.0115 

0.0000 

.0164 

.0346 


-.0266 

-.0227 

-.0147 


.0165 

.0382 


ALPHA  = 


-.0344 

-.0268 

-.0142 

0.0000 

.0175 

.0385 


-.0341 

-.0268 

-.0149 

0.0000 

.0188 

.0407 


Table  XEII-A 

CHANGE  IN  MOMENT  COEFFICIENT  WITH  ELEVATOR  DEFLECTION  ANGLE 


/namic  Coefficients  with  Curve  Fits 


k 


4 
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In  Chapter  2 an  aerodynamic  model  of  the  Jindivik  is  developed. 
The  development  relies  upon  the  existence  of  functional  relation- 
ships between  the  orientation  of  the  aircraft  to  the  air  flow  and 
the  aerodynamic  coefficients.  One  method  of  deriving  these  relation- 
ships is  to  curve  fit  the  tabular  data  with  combinations  of  poly- 
nomials, and  this  is  the  approach  taken  in  this  section.  For 
simplicity,  the  lowest  order  polynomial  which  appears  to  represent 
the  overall  shape  of  the  data  is  used,  and  the  polynomial  coeffi- 
cients are  determined  by  least  squares  curve  fits.  A computer 
subroutine,  PLSQ  (Reference  11) , was  used  to  generate  the  polynomial 
coefficients  for  all  second-order  and  higher  polynomials,  as  well 
as  some  first-order  polynomials  (lift  coefficient  with  flaps  up, 
side  force  coefficient,  and  roll  moment  coefficients). 


Lift  Coefficient 

Figure  A-l  depicts  the  lift  coefficient  data  for  both  the 
flaps-up  and  flaps-down  configurations,  along  with  the  straight  line 
curve  fits  of  the  data.  The  algebraic  relationships  are 

flaps-up  CL  = .08476a  + .132  (A-l) 

flaps-down  = .08476a  + .532  . (A-2) 

Drag  Coefficient 

The  coefficient  of  drag  data  is  shown  in  Figure  A-2  along  with 
the  curve  fits  for  both  the  flaps-up  and  flaps-down  configurations. 
The  algebraic  relationship  for  the  flaps-up  configuration  is 

Cp  = .001506a2  + ,00045a  + .09196  (A-3) 
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HLPHfl  (DEGREES) 
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FIG.  A-2  ORRG  COEFFICIENT  VERSUS  ANGLE  OF  ATTACK 


For  the  flaps-down  configuration,  the  expression  is 


CD  = .001653a2  + .00518a  + .16445 


Pitching  Moment  Coefficient 


Figure  A- 3 shows  the  pitching  moment  coefficient  data  for  flaps 


up  and  flaps  down,  along  with  the  curve  fits.  Tb 


01103a  + .0572 


Cm  = -.001848a  + .0453 


The  flaps-down  curve  fit  is  simply  two  straight  lines 


The  variation  in  side  force  coefficient,  Cv,  is  shown  in 


Figure  A-4  for  flaps  up 


01016  + .00792 


Roll  Moment  Coefficient 


The  roll  moment  coefficient,  C^  (CLL  in  FORTRAN  to  distinguish 


between  lift  coefficient  and  roll  coefficient),  varies  linearly  with 


sideslip  angle,  6,  as  well  as  with  lift  coefficient.  The  wind 


tunnel  data  for  this  coefficient  is  given  in  Figures  A-5  through 


FIG.  fi-S  ROLL  MOMENT  COEFFICIENT  VERSUS 
SIDESLIP  ANGLE 


FIG.  fl-7  ROLL  ilOMENT  COEFFICIENT  VERSUS 
SIDESLIP  AN OLE 
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t coefficie:: 


A- 8 for  the  flaps-up  configuration.  The  defining  relationships  for 


this  configuration  are 


a = .1° 
a = 2.2° 
a = 4.3° 
a = 6.4° 


Cx  = -1.48  X 10-38  - 2.355  X IQ-3 


(A-8) 


C*  = -1.501  X 10-3g  - 1.260  X 10~3 
Ct  = -1.877  X 10-3e  - 1.555  X 10~3 
= -2.302  X 10“ 3 B - 1.445  X 10~3 
For  an  arbitrary  value  of  a between  .1°  and  6.4°,  an  interpolation 
is  made  between  appropriate  curves.  For  the  flaps-dcwn  case,  a 
slightly  different  approach  is  taken.  Since  the  slopes  of  the 
straight  line  curve  fits  vary  with  angle  of  attack,  an  expression 
for  these  slopes  is  formulated.  Thus,  an  algorithm  is  obtained 
that  computes  C\  for  any  angle  of  attack  and  sideslip  angle  combina- 
tion. The  equation  for  the  roll  moment  coefficient  is 

Ct  = Kj_(a)  S + .0012  (A-9) 

where 

K;l(c0  = .1365a  + 1.239  (A-10) 


Yaw  Moment  Coefficient 

Figures  A-9  through  A-12  show  the  variation  in  yaw  moment 
coefficient,  Cn,  with  sideslip  angle  at  different  angles  of  attack. 
These  figures  depict  the  flaps-up  configuration:  the  flaps-down 
configuration  will  be  covered  shortly. 

The  yaw  moment  coefficient  is  a highly  nonlinear  parameter 
which  is  clearly  dependent  upon  both  sideslip  angle  and  angle  of 
attack.  Several  different  polynomial  curve  fits  of  the  given  data 
were  attempted.  While  higher  order  polynomials  result  in  closer 
fits  to  the  data,  fifth  order  curve  fits  are  used  since  they 
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adequately  represent  the  overall  shapes  of  the  wind  tunnel  data. 

Since  no  data  is  available  for  sideslip  angles  less  than  -7°  or 
greater  than  +7°,  straight  line  extrapolations  are  made.  The  slopes 
for  these  extrapolations  are  arbitrarily  chosen  to  be  identical  at 
each  angle  of  attack.  Furthermore , the  slopes  of  the  straight  line 
extrapolations  match  the  fifth  order  polynomials  at  the  "transition 
points".  The  transition  points  (denoted  by  subscript  tp)  are  the 
points  where  the  straight  line  extrapolations  and  the  fifth  order 
curve  fits  merge.  The  general  form  of  the  modeled  relationships 
between  3 and  Cn  for  the  flaps-up  case  is 

3 £ Btp  : Cn  = a~3  + b- 

Btp  £ 3 £ &tp  : Cn  = CqB5  + C2 3H  + C 333  + C432  + C5B  + Cg  (A-ll) 

+ 

3 £ Btp  : Cn  = a+3  + b+ 

Table  XIV-A  below  lists  all  of  the  above  parameters  for  each  value 
of  angle  of  attack.  Linear  interpolation  is  made  for  .1°  £ a £ 6.4°. 
For  angles  outside  this  range,  the  appropriate  extreme  curve  is  used 
(e.g.,  for  a < .1°,  the  a = .1°  curve  for  Cn  versus  3 is  used  to 
determine  Cn) . 

The  data  for  yaw  moment  coefficient  with  flaps  down  is  not 
nearly  as  predictable  as  it  is  for  the  flaps-up  configuration.  In 
order  to  obtain  a crude  model  for  this  configuration,  an  average 
value  was  taken  from  the  values  of  Cn  at  a = .33°  and  a = 2.45°. 

Three  straight  lines  were  used  to  roughly  approximate  the  yaw  moment 
coefficient  versus  sideslip  angle.  The  averaged  values  of  Cn  are 


given  in  Figure  A-13,  along  with  the  curve  fits.  The  equations  for 
the  straight  lines  are 


0 < -1.05  : Cn  = 2.25  X 1O_40  + .003525 

-1.05  £015  : Cn  = -7.75  X 1CT40  + .002475  (A-12) 

3 > 5 : Cn  = 2.25  X 1O"40  - .002525 

Elevator  Control  Coefficients 

Figures  A-14  through  A-16  are  representations  of  the  elevator 
control  data  and  the  straight  line  curve  fits  of  that  data.  The 
curve  fits  are  done  by  hand  since  it  is  desired  that  the  curves 
approximate  the  values  around  the  origin  more  accurately  than  those 
at  the  limits  of  the  elevator  deflection.  The  reason  for  this 
weighting  is  that  the  elevator  seldom  approaches  the  deflection 
limits,  but  instead  typically  remains  within  5°  of  the  neutral 
position. 
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SECTION  II 


Additional  Aerodynamic  Data  Provided 

Some  additional  data  was  obtained  from  Reference  1.  Specifi- 
cally, the  following  dimensionless  stability  derivatives  were  given: 

Ci  > C1  , Cno,  Cn  , , and  Cn5  . Unfortunately,  these  derivatives 

are  available  for  the  basic  aircraft  (i.e.,  without  an  ACRS  installed). 
As  a first  approximation,  these  values  are  used  directly  with  no  modi- 
fication. This  is  claimed  to  be  a reasonable  approximation  for  the 
latter  three  derivatives  since,  1)  the  vertical  stabilizer  is  typi- 
cally the  largest  contributor  to  Cnr  (Ref  16 4.95) , and  2)  the  dis- 
tance of  the  ailerons  from  the  fuselage  precludes  any  appreciable 
effects  of  the  fuselage  on  the  aileron  control  derivatives.  No  claim 
is  made  as  to  the  accuracy  of  the  approximation  for  the  first  three 
derivatives,  however.  Table  XV-A  is  a compilation  of  the  stability 
derivatives  along  with  their  definitions. 

Derived  Aerodynamic  Parameters 

Four  parameters  were  calculated  from  aerodynamic  data  and  the 
geometry  of  the  airframe:  Cl.  , Cjyj^,  Cj_^,  and  C^.  These  stability 
derivatives  are  expected  to  be  significantly  affected  by  the  presence 
of  the  trunk,  since  the  trunk  lies  generally  in  front  of  the  hori- 
zontal tail.  The  horizontal  tail  plays  a major  role  in  determining 
the  values  of  these  four  stability  derivatives,  as  will  be  seen 
below,  and  hence  the  trunk's  existence  must  be  accounted  for. 

Consider  first  the  a stability  derivatives.  An  equation  for 
Cl^  and  C[.j&  (Reference  16)  is 


Table  XV-A 

SMC  STABILITY  DERIVATIVES  FROM  REFERENCE 
(Note:  Units  of  Vrp  are  ft/sec) 


Stability 

Derivative 

Value  or 

Algebraic  Relation 

Definition 

% 

-.370 

p _ 4 9L 

PSbVT  9P 

Clr 

-,06Cj 

+ .459 

p _ 4 9L 

°lr  pSbVip  9R 

V 

-.193 

1 ( VtV 

t 

r - ? 9L 

1 \B923j 

l5A  pSbVj/  9 5a 

s 

- . 125Cl  + .0209 

P _ 4 9N 

PSbVT  9P 

-.15^5 

P _ 4 3N 

^r  PSbVT  9R 

Cn 

^6A 

.0166  - .009CL 

1 

1 

n _ 2 9N 

^6 A ^SbVr?  96a 

(A-13) 


CLS  - VH  £ 

% * 2!;Laa"H  % 5 U-W 

The  constants  are  defined  as 

Ct  - lift  curve  slone  of  the  horizontal  tail 
aH 

n^  - efficiency  factor  of  the  horizontal  tail 
Vjj  - horizontal  tail  volume  coefficient 

xH  - distance  between  the  center  of  gravity  and  the  hori- 
zontal tail  aerodynamic  center 
c - the  mean  aerodynamic  chord 

- the  rate  of  change  of  downwash  angle  with  respect  to 
angle  of  attack 

Detailed  information  on  each  of  these  quantities  can  be  found  in 
Reference  1 6. 

By  techniques  given  in  Reference  1 6 each  of  the  quantities  in 
equations  A-13  and  A-14  can  be  evaluated.  The  results  are 

CLi  = .02465/° 

= -.0604/° 

The  remaining  two  stability  derivatives  can  also  be  evaluated  by 
techniques  described  in  Chapter  4 of  Reference  16.  The  necessary 
equations  are 

°Lq  = 2CLaH  nH  ?H 
CMq  = ~2CLaH  nH  VH 
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Appendix  B 


Salient  Features 
of  the 

EASY  Aircraft  Modeling  and  Analysis  Program 

A computer  software  package,  EASY,  is  used  throughout  this 
report  to  facilitate  programming  of  the  mathematical  models  as  well 
as  to  aid  in  the  analysis  of  those  models.  This  appendix  discusses 
two  main  modes  of  operation  of  EASY,  the  modeling,  mode  and  the 
analysis  mode,  after  which  a source  listing  of  the  model  used  for 
simulations  of  the  Jindivik  slideout  is  presented.  A detailed 
description  of  the  use  of  EASY  is  given  in  the  user’s  handbook, 
Reference  6. 
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Table  I-B 


DESCRIPTIONS  OP  SOTE  OP  TIE  BLOCKS  IN  TIE  EASY  MODELING  PROGRAM 


Block 


Description 


Computes  linear  and  angular  velocities,  linear  and 
SD  angular  accelerations,  and  altitude  by  integrating 

the  six  equations  of  motion. 

Computes  the  aerodynamic  variables  (e.g.,  a,  B , 
p , etc . ) 

Computes  longitudinal  forces  and  moment 
LO  and  Mb)  from  aerodynamic  variables,  gravitational 

force,  and  external  forces. 

Computes  lateral-direction  force  and  moments  (Yb, 
LD  Lb,  and  Nb)  from  aerodynamic  variables,  gravita- 

tional force,  and  external  forces. 

IT  Integrates  the  input  variable. 

Multiplies  input  variables  by  constants  and  sums 
the  products . 


LG  First  order  lag  transfer  function. 

TF  Second-order  denominator  transfer  function. 

Analytic  function  generator  (generates  exponential 
^ or  polynomial  functions  of  time). 

Computes  forces  and  torques  produced  by  engine 
Lh  thrust. 
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Table  II-3 

DESCRIPTIONS  OF  SOME  OF  THE  EASY  ANALYSIS  COMMANDS 


EASY 

COMMAND 

DESCRIPTION 

SEdLATE 

This  command  initiates  a computer  simulation  from 
specified  initial  conditions  to  a specified  final 
time. 

STEADY 

STATE 

This  command  initiates  the  calculation  of  the 
system  steady  state. 

LLNEAR 

ANALYSIS 

This  command  initiates  the  calculation  of  a 
linear  approximation  to  the  model  at  a spec_fied 
operating  point. 

TRANSFER 

FUNCTION 

This  command  initiates  the  calculation  of  the 
frequency  response  function  between  any  two 
specified  points  in  the  model. 

4 


The  Modeling  Mode 


A mathematical  model  of  1)  the  aerodynamic  and  external  forces 
acting  on  a vehicle,  2)  the  equations  of  motion  of  a vehicle,  and 
3)  the  control  system  of  a vehicle,  is  converted  into  a single 
FORTRAN  program  by  EASY  in  the  modeling  mode  of  operation.  A user 
of  the  EASY  modeling  program  provides  information  to  the  EASY  precom- 
piler through  simple  control  statements  that  define  which  of  several 
functional  blocks  are  to  be  used  and  how  they  are  to  be  connected. 
Some  of  these  blocks  are  listed  in  Table  I-B  along  with  a description 
of  the  function  of  each  block. 

In  addition  to  the  functional  blocks,  FORTRAN  statements  may  be 
entered  between  these  blocks.  This  is  a powerful  feature  of  EASY 
since  this  increases  the  flexibility  of  the  model  significantly. 

The  computer  listing  towards  the  end  of  this  appendix  demonstrates 
this  feature  since  many  FORTRAN  statements  are  needed  to  generate 
the  Jindivik  slideout  model. 

The  Analysis  Mode 

The  final  product  of  the  modeling  program  discussed  above  is  a 
binary  object  deck  which  is  in  a suitable  form  for  the  analysis 
program.  In  addition  to  simulations,  several  other  means  of 
analyzing  the  system  described  by  the  modeling  program  are  available. 
Table  II-B  provides  a list  of  some  of  the  analysis  commands  which 
are  available. 

The  Jindivik  Slideout  Model 

A computer  listing  of  the  input  commands  to  the  EASY  modeling 
program  is  presented  in  the  following  pages. 
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